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ABSTRACT 
The stochastic model of the precipitation regime developed by 
Verschuren and Todorovic can be applied to particular climatological 
eee or geographic areas if two parameters can be estimated from 


past records. These parameters, \, and Ags describe the temporal 


1 
distribution of the number of storms and of the total amount of pre- 
cipitation respectively. The present thesis represents an attempt to 
evaluate these parameters and their variability through the year reas 
thirteen climatological stations in Alberta. 

Eleven stations were selected to represent two east-west tra- 
verses across Alberta at the latitudes of Lethbridge and Edmonton 
approximately, and two further stations in southern Alberta were in- 
cluded to evaluate the effects of elevation. While first-order sta- 
tions were desired, their sparse distribution necessitated the use of 
eight ordinary climatological stations. As these latter stations 
report on a daily basis only, the study was restricted to daily time 
intervals. From twenty-five to thirty-one years of record were avail- 
able at each station. 

The theoretical development of Verschuren and Todorovic is 
summarized briefly with emphasis on the basic assumptions. Estimating 
equations for the parameters Ay and Ay are derived and the variability 
of the parameters (characteristic curves) is established by obtaining 
estimates of the two parameters for 71 five-day periods through the 
year. | 

Comparison ee lche characteristic curves for the thirteen sta- 


tions indicates that significant differences in the distribution of 


1ii 











, _ a oY. _ ee Y - 
' 7 fo ry : ro ( + f 
| lene tr ona aes 
7 a» , TORR evetib : : 7 
7 : 3 A hl oh Lapel i re ee ig, ao fe ROE, Py arg 


me : ‘ q &! ] oa tw i ews ha 2a i ole iL 
‘ sve 


ed ‘Heqolaveb ou a 3 aotredigi 2%, ana to Leben at teadsosa, aT 


\ isatgolosemiis walysitseq od babiqqe sd oep at vorvoboT boa aaties 


. 4 = 
‘ —~ - A : Ps 
gozi betaaiszas ed nas atojSseis%hd owl 2 sacta aidgaieien 3 
" 
1 
5 ; 
- --? om «4 = 5 - ~~ ? 
SA IOGaASS 2 S24TIABoD -h ben J meh: BESS it we she ab 
: a = = 4 
2 , =. ; 
-“S7y 70 Jrvcams Letos of2 io & 6je 39 290807 91 (Ol Sar 
¥ af 
—— M y a . : 
i Oo} 39msJ3J3e OA F 1 : Saea.o 7" : WV hi 23 43) 
~ ke 7 
' + © - - ia 
waov aod dj : y 3 ig ez Te fz inva 


a eT | 4 ‘0 @ 


d ¢ Bi 2 Ih ula is . } . oh rs rl ts waa? tbds 


* +4 ’ i | ~ » tay 
- 32a 9% 3 - S 3178 f7¢ Sic vevelg Pied 
7 a oo 
“ ot... 5 c 4 . ° : * a f ‘ 
nmoinombaA iva : wise! to eehbesigal oda te ottetlA seotoe peewy 
mm, U : : 
i] > . J = te if ,° 
4% oct tad rt) rete % i a : “; i ae i 
? mw SL SU 4A Te » 24869 See ge Cy yO& (V+St ai ZOTHGS . 
: 
Me 7 a | 
= To .— >t y 4 of tf ~ = , =T a | oe i ot mt nel ~ 
> / " is . S411 «#0 & LS ‘G od Po: aie BSJIBPLEVS ‘OD > initia O 


= ~~ ~ 4 =» ¢ 3 = =f » f 
¢ e : 1 O03) B22 2-9 Suey o72 «6.elisvystal 
- : 
; notsa rioas ya alde 
=—f 
e 
o> oF ty . 5, 4 f r ; 
7 ati me? a" 1% a i 0; 289 GULSVSsx LG 533% iJ Si 1 ’ 





ae _ an — a 7 = . : ~* 7 “ 7 » 
Ria Puseniseg .arvisgqevess olasd ol) 10 cieadqums fitw vy i%sbs bent aan 






@2iiidetasy of2 bas hbevideb ain A ina 










miatberdo @ fedeliduses st (esvimr -aldat<1s3 s67u/5) tm soneaageads ‘ 
: a _ 7 aT ‘nt - : Tons 


abotaeg Yeb-ovtd. it 202 wasseme%es coz ay. 15 op: 
- 7 7 J av 7 7 - a 7 





of precipitation occur. These differences appear to be partly the 
result of geographic location and elevation and Revie due to different 
operating procedures between first-order and ordinary climatological 
stations. The minimum number of years of record required to provide 
reasonable estimates of the parameters A, and ho is also established. 
Finally, the effectiveness of the theoretical equations is 
illustrated by comparing the observed distribution of the number of 
days with measureable precipitation and of the amount of precipitation 
for varying time intervals with computed values based on the estimates 


of AS and A The reasonable agreement obtained for most stations 


2° 
indicates that the theoretical model is a reasonably close representa- 
tion of the precipitation process. The model provides useful estimates 
of the two main variables, which are the number of days with measure- 


able precipitation and the total amount of precipitation in any 


designated time period. 


iv 
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CHAPTER I 


INTRODUCTION 


1.1 Outline of the Project 


The main theme of this project is to study the precipitation 
patterns at various locations in Alberta and compare their distri- 
butions with those distributions that have been theoretically derived © 
for various variables associated with the precipitation phenomenon. 
In general, this study involves: 

1. The development of characteristic parameters for the sta- 
tion or area under consideration. These parameters will indicate 
seasonal changes and therefore, will have to be considered as time- 
dependent. With the aid of these parameters, the distribution func- 
tions for several variables can be calculated. 

2. The comparison of the parameters for an ordinary climato- 
logical station and that of a nearby first-order station. 

3. The comparison of the characteristic parameters for dif- 
ferent geographic locations in Alberta. 

4. The minimum number of years of record to give reasonable 


estimates of the characteristic parameters. 


1.2 Reasons and Aims 
The full potential of the water resource for a specific area or 


region can only be realized when its magnitude is known. From a 
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hydrologist's point of view, often there is insufficient information 
about the flows in many rivers because of the inaccessibility factor 
or the insufficient instrumentation of such rivers. A commonly used 
technique in estimating river flow is that of the precipitation - 
runoff relationship (Chow, 1964; Linsley and Franzini, 1964;. Foster, 
1949) where the unit hydrograph or other systems are employed. In 
hydrograph analysis, the excess precipitation is converted to runoff. 
Before explicit accuracy in determining runoff can be obtained, more 
emphasis should be directed to the understanding of the occurrence of 
precipitation. 

The seasonal variations of precipitation become important when 
one considers the loss of life and property damage that may be possible 
in the event of a flood or the critical effect of a low water supply 
when a hydraulic structure such as a dam is built across a river. The 
reservoir of any hydraulic structure should be designed to accommodate 
a large storage for the high rate of runoff during spring as well as 
provide ample supply to large drawdown demands during the summer when 
precipitation is less. Computations requiring the determination of 
the amounts of precipitation depend on many factors of which many in 
turn are dependent on the season. 

A common practice in designing the safety margin of a dam 
(Chow, 1964; Linsley and Franzini, 1964) is that of transposing the 
worst storm that occurred in a particular area during a certain season 
to a region that is usually described as climatologically similar. if 
such maximization and transposition of storms is to be allowed, the 


difference and variation of the climate in the two regions with respect 








bee Minomecg A .2revit dave ‘lo elemaemresank: ue acaba 
horsetigise xd ont, to toda, 2¢ wold aSvin Bie AE IAM AAE ‘sunsetiaed 
.rera04 Per sidtizast bis Yelentl amaes wood} qidutetietos Seal 
tl. shevolome ons antet2ye “19nty "0 sateen Jig: “ket ova ewe 
items ot bedrayinga 2: 


yiom ,bohrnddo od neo 2404u5 ania iiss i vovvooed a hotiqay pee 


ho SoavayNone sd2 lo warboitenrshau ef9-o) bejostib 4a dlectd Shee 


Jit aiqeaaay 
nedw ined iogm? emiaad oottss iqrasty Yo. enor rueieN (urgease eral 
kteos od “Bm Seno subesd’ yrtoqeny bis sirl- 26 weet aP embperoo 9m 

(lygva Ietew Wool G to Justis is6'Ttxe aad “<6 sera LDL au tn3vo ity mt 
sy ToVvii-e8 Seatos Sted 2h ited 28-902 ott {2b orlusvkyd ® asd 
ssaspleneae oF botigteeh od Sluode 2¢uPsuate Se onder “iit Jo 1 O-vaeeet 

oi Fiew ds glitqe yaltub Povet us 618" fata 2 wiv 103. ope Tose egual s 

nedw cronme oA aid (Ob ehgayos éobunay Jyteiod “saque Sic ad cbs vor 

Ww cottaniwes seb satan bet Lupoz OE Lagi Ue  seaL as ond ateghneny 

ab ent doidW Yo exoyast ynam 46 bissyeb aod sasietsund: 0 2 Stews, ot? 
‘aaaed oa stg? sasbanageb be aeoulll 
gab. ® te asin Xsothe 5S? “aa or Sattzing’ romp: A: . 


a — 


? 


etiqiosty ce_oxe ans un hans tagger 


» 


B, 


a sant te 9c th eiicoet Has ine een a 


4 
ia 
ia 


7 


a 
7 


» 


3 
to the season of the year must be carefully examined in order to avoid 
improper estimation of the runoff produced by the storm. To reduce 
the complications in developing the design flood through the method of 
transposition of storms, a technique adequately describing a flood 
through the distribution classification of. precipitation can be used 
as an alternative. 

Since no event of the precipitation phenomenon can be predicted 
with certainty, only probabilities of the future outcomes can be deter- 
mined on the basis of the present state. Therefore, those distribution 
functions that are dependent on the season of the year will be most 
applicable to give a description of the variables associated with pre- 
cipitation. 

For any particular season of the year, if one considers the 
number of storms or the total amount of precipitation in a specified 
time interval (variables associated with the distribution of precipi- 
tation), it has been shown (Verschuren, 1968; Todorovic, 1967) that 
distribution functions that are theoretically developed and which make 
use of the characteristic parameters for a station or area can be of 
assistance when precipitation is considered in the engineering design 
of hydraulic structures. It has been suggested by various authors 
(Verschuren, 1968; Thom,1968; Cramer and Leadbeater, 1967) that 
theoretically derived distribution functions allow a better under- 
standing of the conditions for which these distributions hold and offer 
a better chance for explanation and correction of differences between 
theory and the distributions obtained from the observed data than do 


distribution functions obtained by experimentally fitting the observed 
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4 
data to well known distribution functions. 

The goal of any weather modification program (U.S. Dept. of the 
Interior, 1970) is to learn how to manage precipitation in water defi- 
cient areas by cloud-seeding and to do it in an efficient, economic 
and socially acceptable manner. Most programs have been developed 
with the aim of increasing precipitation, especially snowfall in 
mountainous areas, by the cloud-seeding of mainly orographic weather 
systems. Before any programs can be accomplished effectively, an 
extensive statistical study must be undertaken to develop test and 
control areas so that precipitation changes can be accurately determ- 
ined. In order that proper evaluation can be made of the precipi- 
tation at the ground surface, several parameters are required to des- 
cribe the distribution of the various variables associated with the 


precipitation phenomenon. 


Pes Importance of Precipitation in Hydrology 


The emphasis on precipitation in hydrology has arisen because 
records can be obtained easily and cheaply. Aa many regions, records 
have been maintained for long periods and over extensive areas. In 
some countries they are the only available record. Precipitation 
depths can be converted logically to stream flows and hydrographs. By 
appropriate methods, precipitation can also be extrapolated rationally 
to an extreme or limiting value. For example, the probable maximum 
precipitation has many uses in flood estimation. Statistical methods 
are readily applicable to precipitation data. 


The greatest need for statistical treatment of precipitation 
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5 
frequencies, intensities and distribution is felt in various economic or 
related studies. The capacities of sewerage and drainage works are 
usually designed on the basis of a certain depth of rainfall to be 
expected during a specified period of time. Flood control works are 
frequently designed according to a similar criteria, particularly if 
loss of human life is not probable and the protective works are pro- 


posed primarily as insurance against property loss. 


1.4 Difficulties 

When this project was initiated, it was desired to study the 
precipitation distribution in mountainous areas and develop character is- 
tic curves that would be representative of the watersheds of some of 
the major rivers in Alberta. There are many stations in the vicinity 
of the mountains that report precipitation but unfortunately many of 
these are forestry lookout towers and ranger stations with records 
generally ten years or less in duration. Also, a majority of these 
stations are functional in the summer months only. 

At least 30 years of continuous records were ithe Lovee. Ss 
study and the stations selected are operated by the Atmospheric Envir- 
ommental Service of Canada. The only first-order stations in the 
mountainous regions that provide any length of continuous records are 
those of Jasper and Banff. These are inside of the mountain range and 
are considered not to be representative of the river basins that are 
located on the lee side of the mountains. The Ae ecraticne providing 
precipitation data in the vicinity of the mountains are classed as 


ordinary climatological stations and are usually operated by volunteer 
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6 
observers. A description of the different types of climatological sta- 
tions will be given in a later chapter. 

While many stations contained the desired 30 years of precipi- 
tation records, the problem was the fact that the records for some sta- 
tions were not continuous. It was common to find that records were 
missing completely for whole months and in some instances the record 
was missing for the entire year. Another important problem was due to 
the manner in which the precipitation was recorded. For example, if 
the observer noted that it rained on day 1 but did not measure the 
amount of rain until day 2, he would enter a 'C' on the form for day 1 
and the precipitation total would be entered for day 2 and thus would 
represent a two-day rainfall, Pf in tact. ete rained on day 2.) since Jt 
rained on day 1, it was entered as a day with rain. However, it is 
impossible to tell from the records whether or not it rained on day 2, 
but as the observer measured the total rain on day 2, it had to be 
entered as a day with rain also. 

The same type of problem occurs when the observer is absent from 
his station for more than one day. When he arrives from his absence 
and has a measureable amount of precipitation in the gauge, this amount 
is entered in the records on the day of his arrival and is counted as a 
day with rain regardless of whether or not it rained at that time. For 
the days in which he was absent, an 'L' was entered into the record, 
but the days were counted as days without rain as the observer was un- 
certain as to which days the precipitation occurred. 

The above mentioned records are difficult to contend with and 


their use in the analytical computations will lead to futher problems. 
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Therefore, the records have to be removed from the computational work 
and this has the effect of reducing the sample size. This situation 
becomes critical for those stations that already have an abundance of 
missing records and the resultant characteristic parameters and curves 
may not be representative of that station. 

Another important problem that is difficult to control is the 
determination of the amount of precipitation lost through evaporation 
at those stations where the time between observations is not less than 
24 hours. This problem is more realistic in southern Alberta as the 
average temperature during the year is slightly higher than the central 


or northern regions. 
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CHAPTER II 


STATISTICAL TREATMENT OF PRECIPITATION 


ml Description of Variables 


As outlined in Chapter I, the purpose of this thesis is to 
develop a series of characteristic curves that would describe the pre- 
cipitation pattern or behavior for a particular station or geographic 
area. In order to appreciate the statistical approach to the study, 
the main variables are introduced at this time and a brief discussion 
as to their importance is provided. 

The main variables used to develop the characteristic parameters 
are: 

1. The number of storms that occurred during the time interval 
(tt), designated by Nee The time interval for this particular pro- 
ject is five days. 

2. The total amount of precipitation that had fallen during the 
same time interval of (t.»t)- This variable is represented by the 
symbol Xi 

The variable Nt will be used to calculate one of the character- 
istic parameters, Ay (the number of storms per some specified time 
period) while X, is the main factor in the derivation of another 
characteristic parameter, ro (1/the amount of precipitation per storm). 


Further in this chapter, a simplified version of the statistical anal- 
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9 
ysis developed by previous authors (Verschuren, 1968; Todorovic, 1967) 
will show how the theoretical distribution function of the random vari- 
able Nee denoted by F(x/t), and the theoretical distribution function 
of the variable Xs designated by FQ), were derived. 

Other variables used in the theoretical development of the sto- 
chastic process are listed as follows: | 

1, The total precipitation during v storms, X,, where v teat wee 
..., denotes the number of storms. 

2. The total precipitation during the vth storn, Zs where 
v = 1,2,....., denotes the number of storms. 

3. A specific amount of precipitation, x. 

dyeeThesprecipatation intensity atytimert, eee where t denotes 
the length of the time period. 

5. The amount of accumulated precipitation during a single 
storm up to time t, Gee 
2.2 Discussion of the Precipitation Process 

Although the science of meteorology has improved in the past 
several decades, the determination of the outcome of a random variable 
such as precipitation is likely to be done with uncertainty. While the 
statistical approach to the prediction of precipitation will not be 
effective on a short-term basis, it does become useful on long-term 
considerations. The long-term behavior of precipitation would be of 
interest to the hydrologist and hydraulic engineer. As the onset of 
any form of precipitation is dependent on time t with a periodicity of 


one year, then probabilistic laws are likely to describe the behavior 


of future precipitation and the statistical tool to illustrate this 
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10 
behaviour is the stochastic process, 

Precipitation is recorded using a depth in inches and is usually 
in the form of daily, monthly or annual values. Other methods of 
describing precipitation (especially rainfall) is through the means of 
hyetographs and mass curves which are obtained from recording rain- 
gauges, an instrument giving total precipitation during any time in- 
terval. The hyetograph, which is a measure of rainfall intensity (,) 
against time is illustrated in Figure 1. It will be the initial con- 
cept in the statistical description of the important variables. Even 
with today's progressive techniques and scientific models, rainfall 
intensity is one variable that is difficult to predict with assurance. 
Thus future predictions will contain uncertainties. This implies that 
future hyetographs could not be predicted with certainty, thus et be- 


comes a random variable for any t > to: 


Precipitation Intensity g 





Time 


Fig. 1. Hyetograph 
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Statistically speaking, 
Wir Ae eta 


is a family of random variables representing a continuous parameter 
stochastic process for any specific units of time selected from the 
time interval (to >). 

While rainfall intensity is very important in different aspects 
of surface and groundwater hydrology, it is a difficult variable to 
analyze because of its variation in both time and space. In order to 
obtain data for intensity measurements, special instrumentation is 
required (remote-controlled tipping-bucket raingauge) of which there 
are only a small number in Alberta. Thus it would be difficult to 
develop a distribution of precipitation that would be representative 
of the mountain regions or the province as a whole (see Figure Be 

It therefore would be more practical to analyze precipitation in 
the form of accumulated amounts up to some time t, which from Figure 1 
would be the integral of the function oe. Then, Xie the accumulated 
precipitation in the interval (t.»t) can be calculated from the follow- 


ing equation: 
= t 


With xX being a random variable for every. t > t, and because gt is 


{ 5 = 


represents a non-decreasing continuous parameter stochastic process. 
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12 
With a simplified hyetograph as shown in Figure 2 indicating 
the total precipitation during a stormy period to be concentrated at 
the end of that period, one can develop a new function from the process 


scat 


Amount of Precipitation 
During a Single Storm Z 





Time 


Fig. 2. Simplified Hyetograph 


From Figure 2, it is seen that z,, is the total precipitation during the 


Vth storm and this can be derived by taking the integral of Cp» or 


iG 
7 ee OTN eX 
Veta r v 


v-1? 
where X, is the total accumulated precipitation until the end of the 
vth storm. From this process, x becomes a random variable with the 


expression 
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{X, 5 Wien] oes eect 


denoting a discrete parameter stochastic process. 

By transposing the simplified hyetograph into a mass curve as 
indicated by Figure 3, one can establish the number of rainy days or 
periods in a certain time interval. With reference to the stochastic 
process 1X3 t Pape along with the counting of the number of steps on 
the mass curve diagram of Figure 3, it can be said that the number of 


storms (n,) during the time interval (t,»t) is a random variable where 
Ae a 


represents a discrete parameter stochastic process. 


X 
Vv 


Accumulated precipitation 





Time 


Fig. 3. Simplified Mass Curve of Precipitatiom 
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14 
In this study, precipitation will be treated as a purely random 
phenomenon, subject only to periodicity of one year or fractions of a 
year due to seasonal variations. The concepts "days with measureable 
precipitation", "days with precipitation", "rainy days, ''stormy days" 
that are used in this study conform with common usage of these 


expressions. 


2.3 Theoretical Development 


Extensive studies have been carried out by various authors 
(Verschuren, 1968; Todorovic, 1967) in regards to the stochastic pro- 
cess {X,3 te(t,,t)} and reference should be made to those publications 
for a more detailed mathematical treatment. It is only intended in 
this section to show how the theoretical development of the precipi- 
recon phenomenon was achieved by these authors so as to provide a 
guide when observed and theoretical values are compared. 

The important concept of the precipitation variable is that of 
predicting or estimating the average number of stormy days in a time 
interval (t,t). It would be the ability of giving a first approxima- 
tion or estimation of the number of days with rain or snow to be 
expected in a time length of a week or a month during some part of the 
season. Since this concept will be a function of time and other charac- 
teristics of Ny» one has to use the theory of probabilities. It is 


necessary to calculate the following probabilities: 
Pin, So p(t) 


for every t belonging to (t.»t) and v being an integer greater than or 
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15 
equal to zero in order to estimate the number of storm periods in some 
interval (t,t). 
Cant ; 
By letting EY represent the set of all sample functions of 
the stochastic process {X, 5 te(t.,t)I, that is, gee being the event 


that exactly v storms occurred during the interval (t.»t), then 
z Avge 
pct) = PCE O* ) 


Considering that these events are exhaustive and mutually disjoint, it 


follows that: 
z = pcptot 
Pin, v} P(E, ) 
Before the probability of Eva" 4 can be established, several assumptions 
and conditions must be satisfied and one condition states that the 


t,t+At 
lim wee Pike ) 
sa en test 


At+0 AE — Oo 


This implies that the probability of two events occurring in the time 
interval (t,t+At) is equal to zero as At approaches zero. In general, 
it states that only one storm at a time can occur and this will be 
consistent with the usage of daily precipitation or stormy days. For 
the purpose of the analytical approach used in this thesis, a daily 
recording of precipitation will be considered as though it occurred 
from a single storm. 

Another condition that affects the probabilities of the events 


Hon is given by the expression: 
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t, t+At,,ty,t 
jim * SF; Eo oe (Sly fe is ie 
At+0 At 
t,t+At 


In this condition, P(E 7E oop represents the probability that the 


a 
upper bound of the vth storm will belong to the interval of time 
(t,t+At) conditioned on v-1 storms occurring in the interval (t.»t). 


It can be said that the conditional probability depends on t, At and 


vy and can be expressed in the following manner: 


Cirritye tat. 
P(E, /E,0+ ) = A,(t, At, v-1). | 2 


Based on the assumption that the probability that a storm will occur in 
the interval (t,t+At) does not depend on the number of storms up to 


time t, then the function 2.1 can be reduced to the form: 


te t+AG 


P(E, 


) =A, (t, At) 


The above assumption is not completely valid as there is a certain 

~ dependence on the number of storms up to time t, but in order to devel- 
op the theoretical equations, independence was assumed (Verschuren, 
1968; Todorovic, 1967) and this will be consistent throughout the 
remainder of this section. With emphasis on this assumption then, in 
the situation where At is very small, \, can be considered to be a 


1 


linear function with respect to At and the expression 
A, (t, At) = A, (t)At 


is valid. Based on the above conditions and assumptions and with the 


use of generating functions (Verschuren, 1968), it follows that: 
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t 
fA, @) ar] [ 2A (nae 1” 


tost, _ 
P(EO*") = p(t) =e 
vy! 
In other words, 
- f°, (nd ar] [ fay ar 1” 
P(n, = Vv) =e O O s 
v! 
Therefore, if one wished to estimate the average number of storms dur- 
ing the time interval (t,»t), it would be accomplished by taking the 
expected value of Nye or, 
(8 8 Vv 
=¢t A, (x) dr [ te d, (x) dr ] 


E(n,) = vegve | 
v! 


and can be reduced to 
E(n,) = iN (r) dr ae 
t to 1 : 


It is an established fact that the average number of storms or 
stormy days during some specified time interval is dependent on the 
season, thus it can be concluded that E(n,) is a function of time. 
Before Equation 2.2 can be evaluated, the function A, @) must be com- 
puted. Unfortunately, it is an unknown function and therefore must be 
solved from a non-analytical point of view. Let it be accepted that 
A, (©) is always positive (since precipitation can never be negative) 
and that A, @) is periodic with respect to were ral variations. If it 
is assumed that seasonal variations are small over ane time intervals, 
the function 4, @) can be approximated by a constant i, over these time 


intervals. With acceptance of the above argument, Equation 2.2 becomes: 
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Ein eS) dr 
Oe Mee 
and by letting to equal zero, 
E(n,J =A,t - ooo 


By considering a number of years of precipitation records at any 
station or area and with the averaging of the number of storms during 
a time interval (t. >t), a value for Ayt is obtained. The average num- 
ber of storms per unit of time, represented by Aye during the time 
interval (0,t) is found by dividing A,t by that time interval. Hence, 
the probability that exactly v storms occur during the time interval 


(0,t) follows a Poisson distribution and is given by: 


) 
Go ae j 2.4 





v! 


with the mean and variance being A,t when ty is considered as zero. 


From the definition of the events a it was established that 


P(E,07') = Pin, = v} = p(t) 


thus the corresponding distribution function F(x/t) of the random var- 


iable n, has the following form: 


F(x/t) = Pin, cas = Eg? 


where [x] denotes the greatest integer not greater than x. But from 


Equation 2.4, 





p(t) =e vy te 


v! 
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The individual and cumulative terms of the Poisson distribution can be 
found in standard tables or programmed through the use of the computer. 

With the basic understanding of the theoretical development for 
the distribution of N,» One can now establish distributions for other 
variables since some of the arguments, conditions and assumptions men- 
tioned previously will be applicable in further derivations. 

In considering the distribution of x, a specified amount of 
precipitation, the random variable X., which defines the amount of 
precipitation during y storms and is represented by a discrete parameter 
stochastic process, will be employed in the development of the distri- 
bution. Introduce ee to represent the set of sample functions of 
the process {X15 t > al that has exactly v points in the interval 
(x5>x). Then the set GRO will represent the event that the total 
amount of precipitation will be less than or equal to x - Xo during 
exactly v storms and greater than x - Xo during v+l storms. 

From set theory which states that the intersection of events is 
a null or empty set and that the union of events is a full set for 


; vaanee 
mutually disjoint and exhaustive events, then for the set G°’’, 


saat a) S;0" = impossible event 
and 


(ee) 


x k 
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From the above definition then, 


Kok 


P{X, < x} = jer G; 


) = C(x) 


where C, (x) denotes the distribution function of xX for every v= 1,2.., 
Xx ae : 
and the events ae are disjoint and mutually exhaustive. 
Before the probabilities of Ghee can be computed, several con- 


ditions have to be satisfied with the first suggesting that the 


x,x+Ax 
li veo? (GC ) s 
ee 


Ax>0 Ax 
implying that the probability that more than one storm will contribute 
to an infinitely small amount of precipitation is zero. 


The second condition expressed by 


Recpetsyogh 
WG) eee MR 2.6 
Ax+0 Ax 


has the following explanation; the conditional probability 


WoO ne 
1 


end of the vth storm will be between x and x+Ax and is conditioned by 


P(G (G22) declares that the total amount of precipitation at the 
the total amount of precipitation equal to or less than x during the 

previous y-l storms, It is evident that the conditional probability is 
dependent on x, Ax and v and hence can be represented by some function 


that involves these variables, specifically: 
dy = h(x, Ax,v-1) - 


Let it be assumed that the total amount of precipitation during 


one storm is equal to or less than Ax does not depend on the total 
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amount of precipitation during the previous v-1 storms, Then Equation 


2.6 reduces to: 


X,xtAx 


P(G)” ) = Ay Ox, Ax) 


If it is also assumed that do is a linear function with respect to Ax 


when Ax is very small, then 
A, (x, Ax) = A, Ox) Ax A 


By accepting the two conditions and the assumptions and through a rigor- 
ous mathematical procedure involving generating functions and differen- 
tiation (Verschuren, 1968), one finds that, 


X X Vv 
. ee A, (s)ds [ x, i, (s)ds ] | 


P(g.” )ze 


v! 
In this case A, (s) is an unknown, but if it can be assumed that A, (Ss) = 
Ao» a constant, by using the non-analytical argument as in the case with 
\,@), then Equation 2.7 reduces to the fotiowing: 


5 
p(G*or*) = e ee 


v! 





with Xo taken as zero. 


It can be expected that ho (S) is a function of time as the prob- 


Xoe 


¥ is different during selected periods or 


abidaty of occurrence of G 


seasons of the year and . 
Vv 
C(x) + P{X thexbye on 16 iba] Pes 


or: 


mtiougd eat .cimove I-v auotvety oft gabeith 


J 











CdR Cavoys 


“ o% aaagnet Keiw noatoout tevnides, at nS sucd bemuees oely ne 
~ + giora 4 Lae ve ae a ant 


xA(K) A = Gaye A - ’ i 


syogir x dydbsdy tne enoisqnseee silt vie eacitihnos owt ont yaitqesod Ye a 
-(OTS7T LD baw EO bIamMA grisereneg, gnivilvaL oxlbepcnay: txol tanedTamt +o 
op) abel sie ae Jnetioiien sv) iia 





fy ' aK | ~ 
- ia a 4 | wh a | - 7 : 
- f on £e “+ } Ce) eh 4 ; 2 ipk- i . "1 
rt : oat ool sai," ere . = an 
u& i ®y J U t } y 


= fe) & Sour boadlees ot no 77.2 ei few OR A tts att ams {2y, cA seed atte ae 


bur ; 
fitiw sana SHt ai ea Tiroatiyte 14 ia-tiam-now aif? asic xd ine? antes Bagh 


t 


4 pgnewyl lod “eds One eeuber 14 ne csaupel ait 





ly . ‘y . i : 20 ; Po} oe : - 
biel ea ee uk a ; . aa 
Lae el ? |  — i 


eee 2a i 
siaitis ieiks ix ante x0 mobromi © el C)gh 20? festoagia ad min ca 
at wee ‘ 


= me ahoiry berdalee ainianb: tires =e i Tonio sancti 







id 


5 j 
Gos leemS e SEs eee) 


0 a) 
Through further mathematical interpretation of Equation 2.8, 
E(X,) = V/A 


and 


w/d 


Var (X,\) 9 


With the distribution of x derived, it is now possible to show 
how the distribution for Xo which defines the amount of precipitation 
during the interval (tt), was developed. For first considerations, 


choose an event B, (x) such that, 


< x} ¥ x > 0 


B, = {u); X, 


which defines the event BE) as being those values where the amount of 
- precipitation during the interval (t.»t) is less than or equal to a 
specified amount of precipitation greater than or equal to zero. By 
taking the probability of such an event, P{B, (x)}, a one-dimensional 
distribution function of the stochastic process {X,3 lee t} is obtain- 


ed which will be denoted by F(x) where 
Csr P{B, (x)} = PLw; va eels 


On the basis of the set theory of mutually disjoint and exhaus- 


tive events, the event BL (x) can also be written in the following form: 
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Ba ere tet 
B(x) = LU) [E° () 3,01. 


Thus, 
F.(x) = P(B.@)} = JF PLE eae B, (x)]. 2.9 


With the aid of Bayes theorem and the logical interpretation of sets 


and events, Equation 2.9 becomes: 


‘ @ tone eax 
F, (x) iy ey EY (Ss; je 


With the acceptance of the assumptions that Ay and 5 are constant and 
to,t Xo»X ; ; 

that the events EY and ee are independent of each other, it can 

be proven theoretically that the Poisson distribution represents. the 

probabilities of these events (Verschuren, 1968). This allows the 


distribution function of x to be written in the following form: 


v j 
iy oN Ghee KT a (Aa tt) (Ax) 
Bea OS epimcrmate to L 1 z q 2.10 


-For x = 0, Equation 2.10 reduces to 


eet 
Be (O)2ste0 4 


which allows rapid calculations of the probability of zero precipitation 
during the interval (t,t). Ganaires programs were used to calculate 
the values for the theoretical distributions shown in this chapter and 
one program is presented in Appendix F. 

To calculate the average amount of precipitation during a time 
interval of length t, the expected value of xX. is taken, which in effect 


is the expected value of Equation 2.10. After a lengthy mathematical 
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procedure (Verschuren, 1968), it was shown that, 
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and the variance was equal to 
2ait 
Var(X,) =—=>. 
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CHAPTER III 
PRECIPITATION MEASUREMENT AND THE TYPE OF RECORDS IN ALBERTA 


3.1 Instrumentation 

With the scientific advancement in both fields of hydrology and 
meteorology, there is a need for improved instrumentation as well as 
the necessity of having complete and continuous information about the 
elements related to the fields mentioned above. While significant pro- 
gress (Teweles and Giraytys, 1970) has been made in the development of 
cont inuous recording devices in the past decade, most of the devices 
are confined to urban areas ee aie at major airports where trained 
personnel are responsible for their maintenance. Continuous recording 
devices are also located at experimental and other research stations, 
but it should be noted that the length of record is usually short, i.e., 
- ten years or less. For all other communities and regions that report 
precipitation, a standard raingauge is employed and records are main- 
tained by volunteer personnel. 

The official Canadian raingauge is a small cylinder with a cross- 
section area of ten square inches. The gauge is placed in such a posi- 
tion that it is free from all obstructions which might interfere with 
the catch of rainfall. The gauge is usually located on a level grassy 
plot, and the rim of the gauge is one foot above the surface of the 


ground. The rain is caught in the gauge and then measured to 1/100th 
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of an inch in a glass graduate. 

Freshly fallen snow is measured in inches and tenths as it lies 
on the ground. Observations are made as representative as possible by 
averaging several measurements and by avoiding snow drifts and wind- 
swept bare spots. The depth of water resulting from melting newly 
fallen snow has been taken as one-tenth of the depth or thickness of 
the snowfall. Thus the total precipitation for any given period is 
obtained by adding together the total rainfall and one-tenth of the 
depth of newly fallen snow. Early in the 1960's most principal obser- 
ving stations were equipped with snow gauges to measure the actual 
water content, but the data do not extend over a sufficiently long per- 
iod for inclusion in the total precipitation being reported. 

Because of the limited coverage of available raingauge stations, 
not all precipitation data may be obtained from properly constructed 
and exposed gauges. Not infrequently the most intense portion of a 
storm centers over an area where there are no gauges and the greatest 
- rainfall is not measured, or at best is caught in some vessel that may 
act by chance as a receiver. The value of each such measurement must 
be judged individually by the person utilizing precipitation data that 
is not obtained from a raingauge. In some cases he can obtain adequate 
and satisfactory data whereas in others he may obtain only good guesses. 
The probable error of such data, of course, will be appreciably larger 
than that of data of standard gauges, but this reduced accuracy merely 
diminishes without destroying the usefulness of the results. For the 
purpose of this thesis, it is accepted that all precipitation data for 


the stations studied were from the standard raingauge used in Canada. 
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3.2 Type of Stations © 

The order of importance of any station recording meteorological 
parameters is based on the 24-hour period that defines a climatological 
day, the period of record, and the elements being observed. 

First-order stations are maintained and operated by the Atmos- 
pheric Environment Service of Canada and are staffed with trained per- 
sonnel. These stations are located on airports and in larger urban 
areas and are operated on a 24-hour basis with most of the meteorolo- 
gical parameters or elements being reported. Some of these stations 
are equipped with recording raingauges, nevertheless, many of the 
_ gauges have only been installed in the past decade and records obtained 
thus far will not be suitable to statistical treatment. 

There are 25 first-order stations in Alberta that provide a 
high quality in precipitation data and the locations are shown in. Fig- 
ure 4. From the map, it can be seen that only a few stations are in 
the vicinity of the mountains. For those surface-water hydrologists 
- and engineers interested in watershed studies, a lower standard of 
quality of precipitation data must be obtained from other climatolo- 
gical stations. 

At ordinary climatological stations, because the majority of the 
observers at these stations are co-operative observers, the times for 
taking the observation varies from station to station. They are encour- 
aged to take the observation twice a day as close to 0800 and 1700 LST 
as eg eS These times may vary as mich as one hour on either side 
of the time indicated. Some stations take the observation only once a 


day, in the morning afternoon or evening. Among the ordinary climatolo- 
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28 
gical stations, the order of importance in regard to the quality of 
observation is as follows: 

1, Stations operated by the Canada Department of Agriculture. 
At these stations , many of the meteorological elements are recorded by 
standard gauges as well as by some continuous recording devices. 

2. Research stations which are under the joint co-operation of 
both Federal and Provincial Research Councils. Temperature and precipi- 
tation are recorded at these sites using standard gauges only. 

3. Experimental Project Stations which are contracted by the 
Canada Department of Agriculture and which only report several elements 
using stamlard gauges. 

4. The stations operated by the Alberta Forest Service which 
are located at the various lookout towers and ranger stations. Except 
for a few stations, most of these are operated during the summer months 
only and official records did not originate until 1962. Precipitation 
is measured in a standard gauge at least twice a day at the tower sites 
- and only once a day at the ranger stations. 

5. The stations that located in the various townsites, ranches 
and etc. with the people in charge being voluntary observers with the 
minimum of training. 

It should be noted that the highest quality of observation of 
precipitation data would be from experimental watershed studies. How- 
ever, there are a select few in Alberta and these did not commence 
Lay RC Te the International Hydrological Decade was started in 
1965. Figure 5 is a map of Alberta which illustrates the locations of 


all continuous precipitation recording instruments as well as the year 
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of commencement or the years of operation. 


3.3 Recording of Precipitation Data and Errors Involved 


A historical description (Canada Meteorological Branch, 1970a) 
will now be given outlining the different manners in which the climato- 


logical day was defined through the different years of operation. 


3.3.1 First-order Stations 

In the period from 1878 to May 31, 1924 inclusive there was only 
one observation a day, and it was taken in the morning at 0700 LST. 

The climatological day began on day N following the 0700 LST obser - 
vation, ended on day N+1 at 0700 LST, and was credited to day N+1 for 
all observed elements. 

From June 1, 1924 to December 31, 1932 inclusive if there was 
only one observation a day and it was taken in the morning at 0700 LST, 
the climatological day began on day N following the 0700 LST obser- 
vation and ended on day N+l at 0700 LST and was credited to day N+l for 
“all observed elements. If two observations were taken each day, one in 
the morning at 0700 LST, and one in the ‘evening at 1900 LST, the climat- 
ological day began on day N following the 1900 LST observation and 
ended on day N+l at 1900 LST and was credited to day Ntl. 

From January 1, 1933 to December 31, 1940 inclusive if there was 
only one observation a day and it was taken in the aeaang at 0630 LST, 
the climatological day began on day N following the 0630 LST observation 
and ended on day N+l at 0630 LST and was credited to day N for most 
observed elements. If two observations were taken each day, one in the 


morning at 0630 LST and one in the evening at 1830 LST, the climato- 
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logical day began on day N following the 0630 LST observation and ended 
on day N+l at 0630 LST and was credited to day N for all observed ele- 
ments. 

During the period January 1, 1941 to December 31, 1954 inclusive 
most stations took four ap sev ati oie daily at the fixed times of 0130, 
0730, 1330 and 1930 GMT, with some stations taking one to three obser- 
vations. The climatological day began on day N following the 0730 GMT. 
observation and ended on day N+l at 0730 GMT and was credited to day N 
for all observed elements. 

From January 1, 1955 to May 31, 1957 inclusive most stations took 
four observations at fixed times of 0030, 0630, 1230 and 1830 GMT with 
some eeitrdnt taking one to three observations. The climatological day 
began on day N following the 1230 GMT observation and ended on day N+l 
at 1230 GMT and was credited to day N for most meteorological elements. 
For the next four years until June 30, 1961 inclusive, the above cri- 
teria applied for the climatological day, however, the times of obser- 
“vation were taken a half hour earlier, i.e., the new times were 0000, 
0600, 1200 and 1800 GMT. 

From July 1, 1961 to the present most stations take four obser- 
vations daily at the fixed times of 0000, 0600, 1200 and 1800 GMT with 
some stations taking one to three observations. The climatological day 
begins on day N following the 0600 GMT observation and ends on day N+l 


at 0600 GMT with most observed elements credited to day N. 


ae ges Ordinary Climatological Stations 


From 1878 to December 31, 1932 inclusive: (1) if a single obser- 


vation was taken each morning at 0800 LST, the climatological day began 
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onl 
on day N following the 0800 LST observation and ended on day N+l1 at 
0800 LST and was credited to day N+l for all observed elements. (2) if 
a single observation was taken each afternoon at 1700 LST, the climat- 
ological day began on day N following the 1700 LST observation and end- 
ed on day N+] at 1700 LST and was credited to day N+l. (3) if two 
observations were taken each day, one in the morning at 0800 LST and 
one in the afternoon at 1700 LST, the climatological day began on day 
N following the 1700 LST observation and ended on day N+1 at 1700 LST, 
and was credited to day N+l for all observed elements. 

From January 1, 1933 to the present, the above conditions hold 
except for several changes. In the case of (1), the climatological day 
is credited to day N instead of day N+l, and in (3), the climatological 
day begins on day N following the 0800 LST observation and ends on day 
N+1 at 0800 LST, and is credited to day N for most observed elements. 

Because of the changes in the definition of a climatological 
day during the time since the first official record, an error will have 
- been introduced, However it will remain fixed and should not cause any 
concern in the analytical computations. An example of such an error is 
as follows: 

1. At ordinary climatological stations prior to 1932, if it 
rained after 0800 LST on June 30th, the amount of rainfall was read at 
0800 LST on July lst and the total rainfall was registered for July lst. 

2. Since 1932, if it rained after 0800 LST on June 30th, the 
_ guage was read at 0800 LST on July lst. Howevers the total amount was 
credited to June 30th. 


Thus, in the case of this thesis where precipitation is studied in five 

























ae Te oes, 


. Soe ae 

$e Deu, ab fic atin bie iweiryseto “Se Opso tld a 
ti(8) .atonmete bovrerdo Lia tension eli 02 ole 
-taalio’sis .7e COTL te nopertotia signs aertey exe nos tmwreede elyaten - 
he bre woitevrseds 121 VOT! was gaiwol let. ¥ (eb nd nngenigt: Laotgeto 


ows 2k (2) sl4Koyeb. oF bes ebets fan bie Tet GOTT te Crh Nb: dom | a4 





bone Yad GUB0 +2 gnicxom ofy ai ono jyeb-ulone odes oTav enoiaviende 
yah 96 naged Yeh lasigolotamsia ond (Tal OOVs “oa norm ae oad mb ene a g 
TA) .00FL du 14h yah ne bybne Len awidavrsedo Ted OONE: ‘pil priwonteien i 
atriemols bevrsedo Ila Gt Lt " yeh oF ‘ned ihots aew esas 

biGd enoigihnos syids anc: | timesTO eto re veh i) wtaunet test ; 
yet fasitgniogeasio ait ,{i) 2 ial ails wt zegrnid Laoiye) aot Sqene 
tadhgctasamils a63 ft) al bos Tho WA) 3.0 Bags etd rr eli og bet tbena. ed = 
yah na zbhe bee noizivreedo Te) 0980 “ena yniwol lope val “the eniged senb 
sginemeia Laytarto: s2e7 ier Maoh ry, hegibers °2i ban (Te. nee che few ea 
bie ; t= : 
fu abgeeaamito B LO nottrastob tH) at egress, ott 2a szimgad nl Gye Say 
oven Elew reo ne ( OroneT. feijitia. tenet eis sankey add’ ant otra age ; 
(it wEnGd ron Hiuvede bra bexft, siamecTPiW caaeaden. becitborsnd ieee 
2: SON 9ne Ave Fo -ebtihexe 1 2110 Soave igntey ini Pers ee) : - 
| i ae 

oh Rb REO! ot toi enditete inoigonaiumica agate Sih i ake 


Br heer thw ‘Riekaker Yo. - tame. itd ggOe orl, ito Tal oven 1a aes 


age Yet at Lubenaiih eb ioataat bes coat aah a 
mt ms any oo per: Sa hie _ 





= 


32 
day intervals, if the years of record overlap the years in which changes 
were made, then there will be an incorrect count of the number of days 
with precipitation in some five-day intervals. This will not change 
the characteristic curves appreciably, but it will have a slight effect 
on the values for individual months. 

As in other types of obervations, those made of precipitation 
are subject to errors of various kinds and from various sources. There 
are first the accidental errors in measurement of catch. These errors, 
however, are small when compared with the daily variations of precipi- 
tation and may usually be neglected. A more serious type of error is 
that due to faulty observations, a type that is sometimes designated as 
blunders, Faulty observations may involve errors of appreciable size 
and in storms of relatively large intensity may vitiate important con- 
clusions. 

Constant errors are more serious than those mentioned above 
because they are accumulative. They arise from a variety of sources, 

- such as inaccurate measuring devices, improper exposure, or other 
continuing causes. Proper exposure and location of the raingauge with 
respect to adjacent objects are essential for good records of precipi- 
tation. If the gauge is improperly exposed the catch may be augmented 
by water blown off adjacent trees, Pheaenenie or nearby buildings, or 
by splashing from surfaces too close to the receiver. On the other 
hand, overhanging trees or sheltering walls and buildings may prevent 
the appropriate quantity of rain from entering the receiver. The gauge 
should be located away from such aijerce so that the distance at least 


equals the height of the object. 
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The principal biased errors in measuring precipitation are due 
to three cause: (1) loss by wetting the receiver when transferring pre- 
cipitation to the measuring cylinder; (2) loss by evaporation; (3) wind 
effects. 

The first error for most precipitation gauges is in the order of 
.01 inches or less. It may be regarded as an instrumental error and 
can usually be neglected. Errors due to evaporation depend on a number 
of factors; geographical and meteorological, location of the station, 
air temperature, wind, humidity deficit, and instrumental error (design, 
construction material, etc.). According to Russian investigations 
(Wiesner, 1970) carried out for a number of precipitation gauge enene 
the mean error due to evaporation is 3 to 5 percent of the annual amount 
for all kinds of precipitation, and .02 inches or less for an individual 
measurement. This becomes significant in the case of ordinary climato- 
logical stations where the times between observation are at least 24 
hours. 

The maximum errors in measuring precipitation are due to the 
effect of wind. In the case of rainfall, errors depend on the size of 
falling drops and the wind velocity. In all rainfall measurements, the 
effect would have a positive sign, that is, they would lead to a defi- 
ciency of rain in the instrument. In the case of solid precipitation 
(mainly snow), blowing out is observed with an increase in wind speed, 
and blowing in when the wind speed reaches some critical value and 
snow dei eeang takes place. Thus, the errors due to wind are different 
not only in value but also in sign. 


These may be remedied to some extent by the use of shields which 
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34 
reduce turbulence over the gauges. Without shields, errors of 20 per- 
cent have occurred at wind speeds of 10 knots and 50 percent at wind 
speeds of 40 knots for rain, while for snow, pune of 40 percent and 


70 percent have been known at these velocities. However, these are 


adverse conditions and exposures. 
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CHAPTER IV 


ANALYSIS OF OBSERVED DATA 


Awl Length of Record of Precipitation Stations 


Of the thirteen stations selected for analysis in the Province 
of Alberta, seven were chosen along the southern portion of the prov- 
ince to represent the region from the plains to the mountains, the 
distance being roughly 200 miles. Four stations spanning the same dis- 
tance were selected to represent a portion of the central section of 
Alberta, while two stations in the extreme southern section of the 
province will be used to study the effects of elevation differences. 

The choice of the stations was on the basis that at least 30 
years of continuous records were available. It was soon discovered 
that precipitation records for some stations were missing completely 
“causing the sample size for some months to be as low as 25 years. 
Nevertheless, it will be shown in Chapter V that reasonable statistical 
measures can be made with a smaller sample size. 

In Tables I and II (Appendix G), the name, location, elevation, 
period of record and other meteorological information about the el- 
ements observed are given for each station (Canada Meteorological 
Branch, 1967). From Table II, one finds that at a number of stations, 
snowfall makes a significant contribution to the annual precipitation. 


As was indicated in earlier chapters, the highest percentage error is 
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36 
associated with snowfall because of the various factors that affect it 
at the surface. Thus, it is vital that accurate snowfall measurements 
be made. More first-order stations would provide higher quality obser- 


vations, especially in mountainous regions. 


4.2 Discussion of Precipitation Data and Associated Problems 


In the hydrologist's point of view, all forms of precipitation 
data are important, especially data from short duration but high inten- 
sity storms because of their pronounced effect on hydraulic designs. 
The variations of precipitation with time and space complicate the 
matter and demand a diversity of techniques of analysis. One of the 
problems that arises from analyses of precipitation data is that of 
sampling. For example: 

1. with respect to time, sampling is inferior for those sta- 
tions that are equipped with standard gauges or non-recording devices. 

2. with reference to space, sampling is poor for any area 
except in high-density gauge networks associated with experimental 
plots. 

Because of the scarcity of continuous precipitation records of 
any length, daily values were employed in this study. This will intro- 
duce systematic errors rather than random errors in the computations 
where the observed data are used to determine the statistical param- 
eters. These systematic errors are discussed in the next few para- 
graphs. 

Depending on the definition of the climatological day, if a 


weather system caused rainfall to start at 0800 LST on day N and the 

















} 4 - 
. - oun 
7 


f3 ; a - = 
: a] 
: . eed Vite 
* > conan 
yi 955 * 4) ; tel 4 anos ont HOT “a Vs oi) ; 70 Seue sod I reiwedin: wit be ral ee a 
i i : 7 oe 
vondmemietem Listwote otmtiase, tails fort vei 3 ent tata ot 2: 


5 


ail voeteup ‘vetair -sbt vere bidbw gael iets wmbto-seTk? som 8 on ned 
‘ 7 ; = —_ 


.enolge? eeoubedmion. it vital yes “ened: By 


5 
‘ Oo (oct Faere eT 1 AOS 2e 
es ae tal = psec 4 tees al tc gh 
7 ; ey 
7 i 1 ‘ } n* . iy ni 


i. ; 
; ( sis = % 7 atah vila;  . Uree +404q neg ome sab 


5 


; Bw Lites ! ' [ hi nOTG v! ie —i9 sis er Bi ne Sales 
, “yr eee 


) ADH Dre Cet toes r Jimi settev Lia 


TW Ot: be 1) YiseteVib.d biremes tet 
' 
y ‘ 
; d m 
e i Arise CXK a ; i 
= i 
aititerce Toa 
' 
I e ’ ’ Lao wW 
’ is J > Tar" D ie Wp SLE 
's ; 
\ e *, . £ oT a tts ~~ ea 
: t 
, wun Sues v7 beneb-derit 
2 ‘ 7 
e 
i sar fi ‘ S77 30 Viltastece 6) i. OPNeS 


»Origt DLiw eidT fi0l3 iP Lavaltuats Stow 2 ale ‘iieb ,agwel yaw 
So x ; i : "y 

; i ; i oo 
wulteypenon Gdduri ero1ss ingt. fiat? taltat 2tox3s oi teevave es 
Tae i 7 We 

“mane ifoiwieete vy solmroyob 01 hore ets anh Levra edo one - 
_ ? 2 . f 






- 7 TAY ere Slt i, bare toeib ote 21077 aiameneee ona im, 
o ay a ! & 


i? 

7 7 ce: » a we : 
as er. ae 
: ~4 re 02 


y e a? - a - on 
TL Kab Copigoio tami | 
’ - “er >» 


» Ki) 
ects 





7 


37 

rain continued until 0800 LST of day N+2, then 2 rainy days are counted 
even though the original weather system should be classed as a single 
storm, During the summer months when convective storms are most common, 
a situation arises when several showers or thundershowers pass over a 
station during separated time intervals of a day and provide a quantity 
of precipitation. If this occurred during the period of a climatolo- 
gical day, only one reading is taken and it is credited to one rainy 
day even though two or more storms may have occurred. This raises the 
problem of classifying or developing a definition for storms as well as 
establishing a time interval between rainless periods that would decide 
on the number of storms that actually occurred. This becomes important 
with respect to some of the random variables used in the mathematical 
derivation of the theoretical distributions developed in an earlier 
chapter. Such variables would be the time of the end of each storm 
and the number of storms during a specified period. 

For example, the assumption was made that the probability that 
-a storm will occur in the interval (t, t+tAt) does not depend on the 
number of storms up to time t and that the total amount of precipitation 
during one storm is equal to or less than Ax is independent of the 
total amount of precipitation during the previous storms. In this 
study, a day with measureable precipitation is the smallest time inter- 
val used. Thus, if it rained on day 1 and continued into day 2, then 
the probability of rain on day 2 is not completely independent of the 
eeeirisade of rain on day 1. Likewise, the amount of rain on day 2 
is foe independent of the quantity of rain on day 1. In other words, 


the assumptions made in Chapter II are not satisfied. However, these 
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carry-over effects will not extend back for any length of time. 

The ability of classifying storms as well as determining the 
time intervals between storms would be an ideal area for further re- 
search as it would involve the meteorological aspects of precipitation 
processes and storm movements. 

By Canadian standards, a day with measureable rain is one on 
which 1/100th of an inch or more has fallen, and a day with snow is one 
with at least one tenth of an inch of newly fallen snow. Any amounts 
less than the above quantities are registered as traces. A trace of 
precipitation has no appreciable affect on the total amount of precipi- 
tation but does have some impact when it comes to counting the number 
of storms during a specified time period. By climatological standards, 
a day with a trace is considered as a day without precipitation and 
this is accepted in the computational work of the thesis. But in 
actuality, it does mean the under-estimation of the number of days with 
precipitation. 

With the scientific advancement in meteorology, many scientists 
including hydrologists are taking interest in the meteorological pro- 
cesses that cause precipitation. However, the present recording of 
data does not give any information about the meteorological systems 
that caused the precipitation. Precipitation caused by cyclones can be 
expected to come from a different statistical population than for exam- 
ple, precipitation caused by convection or topography. More generally, 
high intensity rains, although they may fit the same type of distribu- 
tion as low intensity rains, probably require different statistical 


parameters. Yet these different types of precipitation are intermixed 
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to provide the statistical parameters. In this study, because the pre- 
cipitation records do not provide this information, no difference is 


made between the different types of precipitation. 


4.3 Discussion of Characteristic Curves. 

Thus far, the assumption was made that A, (t) and hy (x) were con- 
stant in order to develop theoretical distribution Pane Crone ator some 
of the random variables involved in this study. The ability to apply 
the distribution functions as a test of theoretical values versus ob- 
served readings will require that the above assumptions be valid 
throughout the remainder of this project. However, it should be noted 
that dy [represent ing h, (t)] and ho [denoting A, (x)] will be subject 
to seasonal variations in which the periodicity is one year. Since the 
distribution functions for the random variables do not consider the 
dependency on time of the parameters Ay and ho» their values were com- 
puted independently at different times of the year in order to obtain 
the variation of Ay and rs with the time of the year, i.e., the season, 
and the resulting characteristic AY and do curves. 

For any station or area under consideration, the characteristic 
I curve can be computed by using any function or distribution that con- 
tains only the Ay paramter. For example, from Equation 2.4, 

-A,t (a,t)” 


p, (t) =e 
aw 





which defines the distribution of nN,» we can obtain the expected value 


of Nye denoted by 


E(n,) = At 
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and solving for Aye we get 





From Equation 4.1, it can be seen that observed data over any 
arbitrary time interval t can be used to calculate the Ay parameter. 
It is also obvious that as t approaches zero, Ay becomes an instantan- 
eous value. With the use of daily precipitation data in this study, 

t can not be less than one day. While the significance of dy as an 
instantaneous value decreases with a longer interval of t, it has been 
observed from the records ante the longest number of continual days 
with precipitation is usually less than five days. Accordingly, the 
time interval for this study was chosen as five days. Therefore, for 
all those months with 30 and 31 days, ty was selected as beginning at 
0 hours on the lst, 6th, llth, 16th, 21st and 26th while for February, 
which may have 28 or 29 days, the last interval began with ty equal to 
0 hours on the 21st. It is apparent that the sample size is further 
‘reduced, at least by one day for all those months containing 31 days 
and either 3 or 4 days in February. Since the characteristic curves 
will represent a series of average values, there will not be any ser- 
jous deviation on an annual basis. However, for the individual month 
of February, the \, and A, values will be affected by the absence of 
the last 3 or 4 days. 


Table III (Appendix H) indicates how the i’ values were computed 


1 
from raw data. Once the aT values are calculated they can be plotted 
against the time of year. The time average value of Ay Q,) can be 


superimposed across the points based on Equation 2.2; 
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E(n,) = f° Ay (ear. 


This equation refers to the mathematical expectation of the random 
’ variable Ne where E(n,) can be considered as the area under the charac- 
teristic dy curve for any selected time interval. All the character- 
istic Ay curves for the stations studied are given by Figures 6 to 18. 

It can now be seen from the variations in the characteristic 
curves that the ht function is not constant with respect to time and 
that the approximate solutions for some of the distributions that were 
obtained previously by assuming h, (t) constant during the interval 
(t,»t) are less accurate at some times of the year than at others. In 
some instances, a maximum value of Ay occurs in June or early July, a 
minimum is reached in late July and August and a secondary maximum 
occurs during September. In other situations, for example Figure 14, 
the variation in Ay is quite large during the period from the end of 
April to mid-October. Thus, the time interval (t.»t) during that per- 
_iod must be taken much shorter than for example during March to achieve 
the same degree of accuracy in determining the frequency of occurrence 
of some of the random variables. A further discussion of the character- 
istic Ay curves will be given in Chapter V which deals with the compa- 
risons between geographic regions and differences among stations. 

The characteristic 5 curve is developed in the same fashion as 
A, except that 5 can not be calculated directly but must be derived 
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from the expression, 
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which assumes that At and Ao are constant over a short interval of time. 
As Ay was calculated using a time period of t equal to 5 days, the: same 
interval must be used to obtain the value for E(X,). The procedure 
used to calculate re is indicated in Table IV (see Appendix H). Com- 
puter programs used to do the computational work will be discussed 
briefly later in this chapter. 

The A, values are plotted against the eine of year and the graphs 
for the stations under study are illustrated by Figures 19 to 31. The 
characteristic do curve represented by the smooth line will be consider- 
ed the average value if it is again assumed that changes in the pre- 
cipitation due to the season may be rapid at any time of the year but 
that seasonal changes do not occur every day, With the five-day time 
interval being short with respect to seasonal variations, the ro values 
can be considered instantaneous values. 

An examination of the characteristic ry curves reveals the 
following properties: 

1. The minimum values of do generally occur during the summer 
months of June and July. 

2. The maximum values are most pronounced during fall and 
winter. 


3. In some cases, A. is nearly constant showing little seasonal 


2 
effect, whereas in the case of Edmonton (Figure 26), there is a pro- 
nounced seasonal distribution. 

4, The variation with respect to time for some stations is 


considerable, i.e., the scatter of the points on some graphs is signif- 


icant. 
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Further discussion of the characteristic hy curves will be given in 
later chapters. 

In summary then, the characteristic Ay curve can be computed 
from the average number of precipitation days during five-day periods 
and that the characteristic do curve can be derived from the average 
values of the amount of precipitation during the same time interval, 
but the value of Ay is required for this computation. With restrictions 
that no data be used outside the five-day time interval, then from the 
computational methods used in this chapter, Ay and Ay can be approx- 


imately considered as instantaneous values. 


4.4 Discussion of Computer Programs 


The meteorological data for the thirteen stations were obtained 
from the Climatological Division of the Atmospheric Environment Service 
of Canada. The records were compiled on magnetic tape using the card 
4 type format (Canada Meteorological Branch, 1970b) in which many 
meteorological elements were contained in each record. Several com- 
puter programs were written to abstract the necessary Ay and do param- 
eters, the first series of programs being developed for the first-order 
stations and another series for the ordinary climatological stations. 

Appendix A gives the program for converting the raw data for 
first-order stations. At these stations, a trace of precipitation is 
recorded as a 'T', but for the purposes of this study must be converted 
to zero. As there is 24-hour coverage of the weather elements at these 
stations, it is unlikely that any records will be missing unless by 


clerical omission. The program for calculating the characteristic Ay 
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and ho curve parameters for the first-order stations is shown in Appen- 
dix B. This program is based on the fact that there are no missing 
data. Thus, by having knowledge of the initial year of data records as 
well as the length of such records, all time parameters including years, 
months and days can be entered as integer unit values. This allows for 
rapid calculation of the Ay and A, values and therefore a minimum use 
of valuable computing time. A sample of the computer output for such 
a first-order program is given by Appendix C. 

Since it was discovered data records were missing completely 
and due to the type of observational program at ordinary climatological 
stations (Chapter I), a special type of computer program was required 
to correct the raw data. Such a program is shown in Appendix D. In 
this program, a trace of precipitation, denoted by 'T', is converted to 
zero while all other alpha characters such as 'L' and 'C' are trans- 
cribed to numeric bogus data. Many conversions of data may be necessary 
resulting in some loss in information. 

With the data converted into useable form, a computer program 
to calculate the characteristic Ay and ry parameters for the ordinary 
climatological station was developed. This program is presented in 
Appendix E. For those stations that are missing records for whole 
months and years, bogus data are inserted into the program so as to 
make the record sequence complete. During the stage of actual computa- 
tions of Ay and Ao» the bogus data are Konoved and the values are based 
on years with true data. 

It is apparent then, that ordinary climatological stations with 


missing records and which employ alpha characters in the observing pro- 
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gram will, after conversion of data, have smaller sample sizes. This 
no doubt causes a greater dependence on first-order stations to be more 
representative of a certain area or region. 

The program given by Appendix F provides a solution for the 
distribution function of the random variable x (Equation 2.10). With 
the values of I and 5 for any station as well as different time inter- 
vals (t = 5 days, 10 days, etc.) entered into this program, the output. 
will consist of a table containing the following items: 

1. The value of x starting at x = 0 and incremented by 0.05 
until the distribution reaches a value of 0.999. 
| 2. The value of the distribution function F(x). 

3. The increments in F(x) for every increment of x. 

The values obtained from this program will be used in Chapter VI where 
a comparison will be made between the observed and theoretical distri- 


butions. 
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CHAPTER V 
COMPARISON OF CHARACTERISTIC CURVES 


5.1 Introduction 


In this chapter, the characteristic A, and ro curves will be 


tt 
analyzed to determine the following effects: 

1. The difference in the distribution of precipitation between 
geographic areas, for example, the difference between Medicine Hat and 
Coleman or, between the southern and central parts of the province. 

2. The discrepancies occurring among first-order and nearby 
ordinary climatological stations. 

3. The significance of topography or elevation difference. 

4, The determination of the number of years of record that would 
give reasonable results. 

The comparison of stations will be affected by location, topog- 
raphy, elevation and the type of observing programs and it is expected 
that these factors will result in the anomalies of some of the charac- 
teristic parameters. While it may not be feasible to compare stations 
of different types in the mountainous regions, it was the intention of 
the author to provide some insight into the differences of the distri- 
bution of precipitation in those regions. The reader should carefully 
examine any conclusions reached by the author in regards to the effects 


of the mountains on precipitation. The reader should also bear in mind 
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the number of stations that are available for study in the mountainous 
regions and representativeness of any station to the geographic area 


surrounding it. 


Se Comparison Between Geographic Areas 


5.2.1 Southern Alberta 

Consider the group of stations that represent an east-west cross- 
section of southern Alberta as given by Figure 32. Except for the 
rolling terrain in the vicinity of Winnifred, there is a gradual change 
in elevation from Medicine Hat to Fort MacLeod, then an abrupt change 
due to the foothills followed by a sharp rise into the mountains just 
west of Lundbreck. In this section, an examination of the character- 
istic Ay and Ay curves will be made from Medicine Hat to Coleman to 
determine the degree of similarity or difference among the entire seq- 


uence of stations. 


5.2.1.1 Characteristic Ay. Curves 

Generally, the characteristic Ay curves for all seven stations 
along the southern profile have approximately the same seasonal var ia- 
tion, that is, near constant values during the winter months, a maximum 
value during late May and June, a minimum through July and August and 


a secondary maximum value in the early part of September. At some sta- 


tions, there is a secondary minimum value in the month of October. 


An inspection of the characteristic dt curves shown in Figures 
6 to 12 reveals that the best similarity between stations for the en- 


tire year are those of Medicine Hat and Lethbridge. There is some 
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similarity among the other stations, especially during the summer months. 
The one station that has a distinct departure from the rest is that of 
Seven Persons. It should be noted that both Medicine Hat and Lethbridge 
are first-order stations while the remaining are classed as ordinary 
climatological stations. It is noticeable that the AY values for the 
two first-order stations are greater than for the other stations, 
especially during the winter months, and this can be attributed to the. 
following reasons: 

1. Most first-order stations are staffed with trained personnel 
and a 24-hour coverage of the weather elements is provided. 

2. Precipitation readings are taken every 6 hours resulting in 
4 observations for every climatological day. This means that more 
storms, especially the convective type, are likely to be recorded than 
at eres climatological stations where only one observation is taken. 

3. The advanced instrumentation at first-order stations will 
allow the recording of storms that contribute only small amounts of 
‘precipitation. 

A comparison among first-order and ordinary climatological sta- 
tions will be made later in this chapter and certain conclusions will 
be drawn at that time. 

To determine the effect of the foothills and the mountains on 
the distribution of the number of days with measureable precipitation, 

a comparison was considered between Coleman and two other stations, 
namely, und break and Fort MacLeod. Since Lundbreck was located ina 
valley between the foothills and the first range of mountains, the 


comparison was made between Coleman and Fort MacLeod as shown in Figure 
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33. This will enable the reader to determine the contrasts between 
areas enclosed in the mountains against those areas along the foothills. 
From Figure 33, it is apparent that there are more days with precipi- 
tation at Coleman during the period of October to January while in the 
period from March to the end of September, Fort MacLeod shows a slightly 
higher value for Aye 

Without discussing the broad field of Reet aarerre circulation 
and storm movements, it can be presumed that the flow pattern during 
the winter months provides more days with precipitation within the 
mountains as shown by the curve for Coleman. In the summer months, the 
flow pattern provides more days with rain or snow on the east side of 


the mountains as indicated by the Ay curve for Fort MacLeod. 


5.2.1.2 Characteristic A» Curves 

A study of the characteristic do curves (representing 1/amount 
of precipitation per storm) for southern Alberta given by Figures 19 
to 25 reveals that for most stations, a maximum occurs during the winter 
OnnT a, a minimum value is reached in the month of June and a secondary 
minimum happens between late August and September. As was noted in the 
comparison of the Ay curves, a degree of similarity exists among the 
Ao curves for the two first-order stations of Medicine Hat and Leth- 
bridge while some similarity is apparent for the remaining ordinary 
Ciinatoloedeal stations especially during the summer and early fall 
months. The larger ry values at the first-order stations can be 


accounted for by the fact that more days with measureable snow or rain 


greater than a trace is possible with continuous observations. The 
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50 
explanations given for the departures in the characteristic Ay curves 
among the different stations applies as well to the hy curves. 

It is interesting to note that for some ordinary climatological 
stations, for instance, Coleman, Lundbreck and Seven Persons (Figures 
20, 24 and 25, respectively), the ho curve is almost constant through- 
out the entire year. This implies that the same amount of bar De ieee - 
whether it be rain or the rainfall equivalent of snow, can be expected | 
during any one stormy day or generally, any period of time. In moun- 
tainous country, especially at higher elevations, snowfall can contri- 
bute to the total annual precipitation as much as rainfall, and with 
reference to Table II, it can be seen that the amount of snowfall at 
Coleman and Lundbreck when converted to rainfall equivalent is slightly 
below the annual rainfall. However, in the case of Seven Persons, the 
annual snowfall is just slightly more than one-half of the annual rain- 
fall. This suggests that the ro curve for this station is not represen- 
tative of the conditions in this area. Also from Table II, the average 
annual snowfall for Medicine Hat (48.7 inches) is nearly the same as 
that of Seven Persons (47.6). However, when comparing the ho curves 
given by Figures 19 and 20, there is a significant difference between 
the two stations. On the basis of the differences in the i curves as 
well as the data given in Table II, a comparison between first-order 
and ordinary climatological stations was necessary and will be given 
ier this chapter. 

Having discussed the effects of the mountains on the number of 
days with measureable precipitation by comparing the dt curves for 


Coleman and Fort MacLeod, it can now be determined what changes might 
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5i 
occur in the amount of precipitation per storm by analyzing the Ao 
curves for the same two stations. With reference to Figure 34, the Ao 
curve for Coleman is nearly constant with an average value of 4.0 where- 
ABbip the case of Fort MacLeod, the Ay values vary from 8.2 in January 
to 2.8 in June. From this graph, it can be assumed that the amount of 
precipitation per storm inside the mountains of southern Alberta (using 
Coleman as a representative station) is almost constant from one month 
to the next showing some independence of the flow pattern. On the lee 
side of the mountains, using the Ay curve for Fort MacLeod, one finds 
that the precipitation per storm during the winter months is constant 
with a value of approximately 6.5. However, towards the latter part of 
May and most of June, the ho value is at a minimum suggesting that the 
maximum amount of precipitation per storm occurs at this time. It is 
often the case where storms moving through southern Alberta or in the 
vicinity cause upslope or orographic conditions to extend from the 
plains region to the first range of mountains. As a result, extensive 
precipitation can be expected with the heaviest amounts in the foothill 
areas west of Lethbridge as indicated by the lower values of ry at Fort 
MacLeod and Lundbreck. With a constant Ao value for Coleman during the 
summer months, it is apparent that the orographic precipitation in this 
part of the country has its greatest variation on the lee side of the 


mountains. 


5.2.2 Central Alberta 
A group of stations as shown by Figure 35 were chosen to repre- 


sent a portion of central Alberta where the distance from the mountains 
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was similar to the stations of southern Alberta. It also was important 
that these stations have at least thirty years of record. It is inten- 
ded in this section to compare the characteristic curves of the central 
Alberta stations and then determine the contrasts between different 
zones of Alberta which will be considered in succeeding sections. 

An examination of Figure 35 reveals that the graph is similar 
to Figure 32 inasmuch that the slope of the land from Edmonton to about 
20 miles west of Edson is gentle, then rises abruptly to form the foot- 
hills and is very steep on the west side of Entrance. It may be noted 
that the first range of mountains in this part of the province is at a 
higher elevation than in the southern sections and that Jasper is 
entrenched in a deeper north-south valley as compared to Coleman's east- 
west valley. As in the discussion for southern Alberta, the character- 
istic Ay and Ao curves will be examined from Edmonton to Jasper in order 
to develop some contrasts as to similarity or reveal any significant 


departures. 


5.2.2.1 Characteristic Ay Curves 

The behavior of the characteristic Ay curves for the four sta- 
tions along the central zone indicates a distinct variation for the 
different seasons of the year. Most pronounced is the maximum value 
that occurs during mid-June, July and in some instances, the early part 
of August. Minimum values appear during the spring months and for 
Edmonton only, in October. At Edson and Entrance, a rather constant 
value prevails during the winter months. 


A comparison of the characteristic aT curves illustrated by Fig- 
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53 
ures 13 to 16 indicate that the best similarity exists between Edson 
and Entrance although Edson has a slightly higher value of Ay dur ing 
June and July. With reference to Figure 35, one can see that the range 
of the foothills lies between Edson and Entrance. During the months of 
June and July, many convective storms originate on these foothills and 
are usually fully grown by the time they reach Edson thus produc ing 
more precipitation at this station than possibly at Entrance. This 
would account for some or if not most of the differences between the AL 
values for June and July for the two stations mentioned above. 

The degree of similarity between Edson and Entrance is of partic- 
ular interest as Edson is a first-order station while Entrance is clas- 
sed as an ordinary climatological station although Edson did not become 
a first-order station until 1960. Considering the distinct departures 
that occurred between the different types of stations for southern 
Alberta, it can only be concluded that: 

1. The area around Entrance has the same slope, exposure and 
orientation to peather storms as Edson. 

2. The volunteer observer at Entrance was conscientious and 
reported all storms producing precipitation. 

3. The characteristic Ay curve for Entrance is not represen- 
tative of the area and should be either higher or lower depending on 
the manner in which weather storms affect this area. 

4. The period of record used for Edson was prior to the station 
becoming a per eee station. It therefore is reasonable to accept 
that any similarity would comply with the fact that both stations were 


ordinary climatological stations before the year 1960. 
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To determine which of the above points is most prominant would 
require further research and therefore will not be pursued in the course 
of this thesis. 

A continued examination of Figures 13 to 16 reveals that both 
Edmonton and Jasper (first-order stations) show a considerable variation 
throughout the entire year. However, there is a good resemblance among 
all four stations during the months of May to September inclusive. 

This implies that on the average, most of the central region, at least 
west of Edmonton, has the same number of days with measureable precipi- 
tation during the summer months. During the winter months, say Noven- 
ber to March inclusive, it appears that more days with snow or rain 
occur inside the mountains and across the plains region as compared to 
the foothills area. 

To get a better understanding of the effects of the mountains on 
the distribution of the number of days with measureable precipitation, 
the Ay curves for Jasper and Entrance were compared against each other 
as demonstrated by Figure 36. From the graph, it can be seen that from 
the period of mid-March to early September the two stations have about 
the same number of days with precipitation. This suggests that the 
mountains have little effect on the weather storms passing through this 
part of the province. During the fall and winter period (mid-September 
to mid-March), there is a marked difference in the Ay values between 
Jasper and Entrance. This indicates that the mountains have a signif- 
icant effect during this period resulting in more days with measureable 
rain or snow at Jasper. As mentioned earlier in this chapter, without 


becoming involved with the mechanics of meteorological circulation and 
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storms, it can be presumed that the dominant flow pattern over this 
particular area during the fall and winter months provides more days 
of meteorological activity within the mountain range as compared to the 
lee side. This is demonstrated by the constant low values of Ay at 


Entrance during the fall and winter season. 


5.2.2.2 Characteristic A» Curves 

Similar to the graphs displaying the characteristic dy curves, 
there is a noticeable seasonal variation in the do curves for those 
stations across central Alberta. Generally speaking, from Piette 26 
to 29, there is a maximum ry value during February and March, a minimum 
value extending across June and July and a secondary maximum dur ing 
October. It can be deduced from these graphs that the bulk of the ann- 
ual precipitation occurs in the summer months in the form of rain while 
lesser amounts are contributed by snowfall. 

Because of the topography of the area between Edson and Entrance, 
it was precluded that there would be differences in the amount of pre- 
pieiration per storm. However, the Ay curves for the two stations 
(Figures 27 and 28) were almost in total agreement with each other 
which is significant since there was a high correlation between the AL 
curves. This in effect suggests that both the AL and ro curves for the 
two stations could be combined and the average curves obtained might be 
representative of an area extending from the lee side of the mountains 
to some location between Edson and Edmonton. 

A further inspection of Figures 26 to 29 reveals that for the 


period of May to September inclusive, there is a good correlation be- 
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tween Edmonton, Edson and Entrance. This implies that any meteorolo- 
gical activity covering the area between the three stations is apt to 
produce the same amount of precipitation for any time interval in the 
summer season. In the winter months, except for Edson and Entrance, 
there is a significant difference among the stations. The ho values 
for Jasper from April to September are slightly higher than for the 
remaining stations and this indicates that smaller amounts of precipi- ~ 
tation are produced per storm within the mountainous areas. To obtain 
a better concept of the effects of the mountains on the ho values, a 
comparison was made between Jasper and Entrance as illustrated by Fig- 
ure 37. 

Referring to the elevation chart of Figure 35, one finds that 
the elevations of Jasper and Entrance are similar. However, a shar p 
and high mountain range separates the two stations. From a meteorolo- 
gist's point of view, this would suggest that those storms approaching 
from a westerly direction would leave most of its precipitation within 
the mountains, whereas those storms causing an easterly upslope flow 
would have most of their precipitation left on the lee side of the 
mountains. An inspection of Figure 37 shows that except for January, 
Entrance receives more precipitation per storm for the entire year. 
While the Ay values for Jasper from mid-September to mid-March were 
much higher (Figure 36) than at Entrance for the same period, that is, 
there were more days with measureable precipitation in the mountains, 
the do ricee for Jasper were also larger. The smaller amounts of pre- 
cipitation at ible suggest the ppg orgs that the mountains in this 


area have little effect on the storms erave ling from a westerly danbes 
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evi 
tion or that Jasper lies in a so-called rainshelter where most of the 
precipitation is deposited either upstream or downstream of the station 
location. Extensive studies of the meteorological circulation and storm 
movements would have to be undertaken to discover the impact of the 


mountain range on the do curve in this particular area. 


5.2.3 Contrasts Between Southern and Central Zones 

There will be situations where hydrologists, engineers and other’ 
environmental scientists will be concerned about the distribution of 
rainfall over a large area or across the entire length of a major river. 
If important diversion of water were necessary to augment supplies to a 
populous region, then the rainfall distribution between specific areas 
may become the prime subject of research. If such situations did arise, 
it would be essential to compare various stations of a specific zone 
and eventually develop a statistical mean value of all stations. This 
mean value of the characteristic parameters would approximate the condi- 
tions of that particular zone. With knowledge of the similarities or 
Greterenees between areas, the engineers in water management could jus- 
tify the need to transport water, if necessary, and determine the extent 
and the amount of water that would have to be transported. 

In this section, the stations of southern and central Alberta 
will be combined into a general set of characteristic curves and a 
discussion will be made between the different sets of curves. The read- 
er should be cautioned that the combination of first-order and ordinary 
climatological stations may result in unrealistic values and the valid- 


ity of the characteristic curves for the stations used will be left to 












" . S 
49 3 dg0omn, ote 1 aint et bo ileo-Gec a nl esti. 7s9e85 +682 20 
‘ tf 7 ; » 
(tot 2d ig ) sto i i sjadieanete i : web. ak aot tuatgh x a 
ig : 


: ; ; ; 
; , ‘tte Be? _ Pr 
rf cays : ' 3a : iia (ts 3 “4 rolseaol 


= =a - = 
f 4-9gf) P 2 y vou b lac yf! aJnstey ve 
} 4 boi 
] y* 4 4 iJ 
ro) 
a 
Pe OA 
f IO i, Lac COecmOT. VRS. 
y ; 5% 
: B , Li alalet 
_ - 
Sew ii Opa ? cat 70qam iI : 


% tuo a sHOLUGOg - 






aii -.a9 oe S aise inoys?Ftb sift avewlad ohea ad Llin molest path 
; 4 , 

= ‘27 ail 

o30 be qb30+40 22 d ae AOL S My Lee hed sis am benat guns os Py ms & 





r§ 7 

seIRY Stayt evar wk sluges yam eononae, inake toa tamk 
Ff Foe peu cs hate hale 

=< 7 / eS 

: Lae «© 5; 7 J “ 7 ; bs , ry MG " > ¥, 

et re. pete - y i. aa 7" ah 1 ee heel 


_ 


58 
his judgement. Also, the combination of stations in the mountainous 
areas with those stations on the plains region will result in misleading 
information. Nevertheless, it is only intended to show whether con- 
trasts do exist between geographic zones and to what extent the informa- 


tion could be used. 


5.2.3.1 Characteristic Ay Curves 

Considering the southern Alberta zone initially, the statistical 
average and standard deviation of the characteristic ht values of all 
seven stations combined is given by Figure 38 and it is a eacely 
noted that significant variations occur from May to November. The 
graph suggests that more days of measureable precipitation take place 
during the summer months indicating that the bulk of the annual precipi- 
tation for this particular zone is in the form of rainfall. An in- 
teresting note about Figure 38 is that during the winter months of 
December to April, the standard deviation is fairly large. This may be 
due to several reasons: 
. 1. The large variations in the i values among the first-order 
and ordinary climatological stations will cause an overall lower average 
Ay curve for all stations combined as well as a significant deviation. 

2, Some snowfall is sublimated in the foothills area because of 
the chinook winds that prevail during the winter months. As a result, 
those storms producing small amounts of snow may not ‘be recorded by the 
stations that report precipitation only once a day. 

3. Partly related to point 2, the lack of proper instrumenta- 


tion for the recording of small snowfall amounts at ordinary climatolo- 
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a9 
gical stations will cause measurement differences among stations in the 
same geographic area. 

4. Combination of the stations in the vicinity of the mountains 
with those of the plains will cause some deviation as the mountain sta- 
tions receive more days with snow because of their higher elevations. 

An examination of Figure 39, which is the combined average and 
first deviation of the characteristic i values for the four stations 
in central Alberta, reveals that generally there is a marked seasonal 
variation for the entire year. Minimum values are reached in the months 
of April and October while the maximum values extend across most of 
June, July and August. This summer maximum may be an indication of the 
high frequency of storms (whether orographic or convective) that happen 
in this portion of the province. Another feature that is noticeable is 
the small deviation of the curve for all stations during the period 
from May to September inclusive. This suggests the following items: 

1. On the average, most orographic storms are large enough in 
areal extent to cover a large portion of the central zone. 

2. The number of convective showers or thundershowers are even- 
ly distributed among these central Alberta stations. 

In the winter months of October to March inclusive, there is a 
large deviation from the average curve. Since the combination of sta- 
tions involved both first-order and ordinary climatological stations, 
it would appear that the deviation in the number of cae of measureable 
precipitation (generally snow) during that winter period is mainly due 
to the type of observing program, topography and storm path movements. 


In order to study the contrast of the number of days with rain 
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60 
or snow for both southern and central Alberta, one must refer to Figure 
40 which displays the characteristic curves for the regions concerned. 
It is plain to see that the combined stations of the central zone 
receive more days with measureable precipitation than the southern part 
of Alberta except for April where the Ay values are similar. The con- 
trast of any significance between the two zones occurs during June and 
July. In that interval of time, the southern region indicates a de- 
creasing value of hy and the central region has its maximum values. 
This contrast may be partially explained by the higher frequency of 
convective showers that take place in the central areas through the 
summer months. Both regions experience a minimum number of days with 
precipitation in October and it appears that the winter storms affect 
the northern half of the province more often than the southern half if 


the hi curves are any indication. 


beeease Characteristic do Curves 

With some knowledge of the variations of the characteristic AV 
Pures and the reasons thereof for the two regions under study in this 
section, one can now examine the characteristic Ao curves for the same 
regions and thus evaluate the quantity of precipitation per storm. 

Considering the southern half first, one notes from Figure 41 
that there is some seasonal variation of the ho curve, but it is not as 
pronounced as the hy curve. Although a minimum yale of approximately 
4.0 exists during June and September, the curve shows constant values 
for the period of May to September inclusive. During the cooler months 


of October to April, much of the precipitation is in the form of snow 
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61 
and higher dy values indicate smaller amounts of precipitation (rainfall 
equivalent of snow) per storm. The large deviations during those cool- 
er months can be explained by the same reasons given for the deviations 
in the Ay curve where the measurement of snowfall, the averaging of sta- 
tions from the mountain and plains areas, and the different types of 
observing programs are considered. 

A study of Figure 42 shows that the seasonal variation of the 
characteristic Ao curve for the central stations is similar to the 
curve given in Figure 41. PheeeotSeane values from June to August 
suggest that, on the average, most of the area from Edmonton to Jasper 
receives the same amount of precipitation per storm whether it be of 
the convective or orographic type. It is interesting to note from Fig- 
ure 39 that the Ay curve for central Alberta is at its maximum during 
the same period and witha little smoothing, a constant Ay value can be 
obtained. This suggests that the AL and ho curves for central Alberta, 
shown by Figures 39 and 42 respectively, would be representative of an 
extensive area west of Edmonton for the three summer months. Another 
feature of Figure 42 is that in the month of February, the least amount 
of precipitation can be Me ebtdatoree the central region. The signif- 
icant deviation over the winter months is due to the lesser amounts of 
precipitation received at Edmonton and Jasper as compared to Entrance 
and Edson where topopraphy has a greater effect. 

It is interesting to note from Figure 40, which gives the compa- 
rison of the Ay curves for southern nd central Alberta, that the dif- 
ferences between the two curves (except for April and May) were quite 


large. But, when the do curves for the two zones (Figure 43) are com- 
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62 
pared, the variation between the two curves is not as large. Using the 
values from Table II, the combined average annual precipitation for the 
southern stations is 16.58 inches and is 18.86 inches for the central 
stations. The variations in the Ay and dy values for the two areas 
would account for the differences in the annual average precipitation. 
A specific study of Figures 40 and 43 shows that the Ay values for 
April are similar, but the ro curves indicate van ener Case ofe2:, Saror 
the same month. Also, in the period of mid-July to mid-August, the Ay 
values are identical but a large departure exists between the aT curves 
for the same period. 

Even though the two areas of central and southern Alberta have, 
on the average, the same number of days with measureable precipitation 
during the month of April, the southern region receives more precipi- 
tation per day. This could be due os the greater amounts of wet snow- 
fall in the mountains and more rainfall across the plains area result- 
ing from slightly warmer temperatures in April over southern Alberta. 
= The other interesting anomaly is the situation in late July and 
early August where the greatest difference occurs between the Ay curves 
yet the Ay values for that period are similar. During this time of the 
year, a larger number of the storms are convective in nature and it 
appears from Figures 40 and 43 that while the central area of the prov- 
ince receives approximately twice as many storms, the amount of pre- 
cipitation per storm is nearly the same. It can only be presumed that 
the fewer Abe across the southern regions provide a greater abundance 
of precipitation or that the type of observing program has a serious 


effect on the characteristic curves of the stations in these regions. 
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63 
Some insight to the above mentioned anomalies may be revealed when a 
comparison is made between a first-order and nearby ordinary climatolo- 
gical station later in this chapter. 

After discussing the characteristic curves for the different 
zones in Alberta using the average waite of a series of stations in 
each zone, it might be of interest to examine the curves for individual 
stations of each zone, especially near mountain areas where the majority 
of the rivers in Alberta have their origin. Referring to the elevation 
charts in Figures 32 and 35, one notices that Coleman and Jasper are 
within the first mountain range (looking from the east) and that the 
stations of Fort MacLeod and Entrance lie on the east side of this 
first mountain range. 

Geroaiening the lee side of the mountain first, one can see by 
Figure 44 that the characteristic Ay curves for Fort Macleod and En- 
trance show PtAEGerenae during the period from early June to mid- 
September and are reasonably similar for the remainder of the year. In 
the early summer months, major storms affect the Province of Alberta 
and usually last for several days providing a significant amount of 
precipitation. From Figure 44, it is demonstrated that Entrance has 
more days with precipitation for the month of June. This suggests that 
those storms affecting Alberta during this time have a tendency to last 
a little longer in the central sections of the foothills. 

In July and August, convective storms are more common and the 
larger ht values at Entrance suggest that a higher frequency of these 
storms occur in the central parts of the foothills area as compared to 


southern sections. The frequency of the convective storms will be 
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- 64 
dependent on the amount of vegetative growth in the area, evaporation 
and transpiration and the weather systems moving through the province 
during those two summer months. For the remainder of the year, it 
appears that both stations record nearly the same number of days with 
measureable precipitation. 

When making note of the amount of precipitation per storm for 
both Entrance and Fort MacLeod, it is surprising to discover from Figure 
45 that the ho curves are nearly identical. This gives rise to several 
suggestions: 

1. The similarity is purely coincidental. 

2. The two stations receive the same amount of precipitation 
per storm. 

3. Because of the departures in the Ay curves, the differences 
of exposure, orientation and topography between the two stations may 
have a significant effect on the do curve, making them similar. 

Discounting point 1 as being by chance and accepting points 2 
dnd 3, one can say that generally on a long term basis, the lee side 
of the mountains receive the same amount of precipitation per storm 
whether it be southern or central Alberta. Point 3 would require 
further research to determine the effects of exposure, orientation and 
topography on the characteristic do curve. 

With some knowledge of the distribution of precipitation on the 
lee side of the mountains, a.study can also be made about the distri- 
bution le the first range of mountains. It should be remembered that 


the effects of topography on the characteristic curves of mountain sta- 
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65 
tions will be important and any comparison made may be unrealistic. 
However, from a practical point of view, the comparison of mountain sta- 
tions may provide some information about the distribution of precipi- 
tation in this difficult area of study. The townsites of Coleman and 
Jasper were selected and the characteristic Ay and do curves for both 
stations are given by Figures 46 and 47, respectively. Coleman is an 
ordinary climatological station while Jasper is classed as a first-order 
station. 

An examination of Figure 46 indicates that except for the period 
of mid-March to the end of May, Jasper has considerably larger Ay 
values than Coleman, but by Figure 47, it also has greater Ay values. 
This suggests that while it has more days with precipitation than Cole- 
man, the amount of precipitation per storm is less. It is obvious from 
Figures 46 and 47 that Coleman receives more precipitation during the 
period of March to May inclusive and this is confirmed by Table V which 
gives the 30-year monthly precipitation for those months mentioned 


above. 


TABLE V 


30-YEAR MEAN MONTHLY PRECIPITATION (INCHES) FOR COLEMAN AND JASPER 


Station Month 

March April May 
Coleman 1.26 1.56 2.02 
Jasper . L100 0.99 0.63 


Several interesting deductions of the distribution of precipi- 


tation inside the mountains can be made with reference to the geographic 
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66 
zones of Jasper and Coleman, the characteristic curves of Figures 46 
and 47, and Table V. Some are as follows: 

1. Higher Ay values at Jasper indicate the passage of more 
storms through that region. 

2. Large values of hy for Jasper suggest small amounts of pre- 
cipitation from the storms indicat ing that the townsite is in a protect- 
ed area (meteorologically speaking). 

oe Fewer storms moving through the southern mountain areas pro- 
vide more precipitation as compared to the central regions. 

4. Coleman, being at a higher elevation will likely receive 
more precipitation and may not be representative of the southern moun- 
tains in general. 

5. The effects caused by the different types of observing pro- 
grams in relation to measurement of precipitation. 

From the factors listed, it is obvious that more stations should 
have been selected within the mountainous areas and greater emphasis 
placed on topography, location, orientation, exposure and meteorolo- 
gical circulation of the atmosphere. Limitations were imposed on the 
number of stations chosen for analysis and an involved meteorological 


aspect was not the scope of this study. 


5.3 Study of First-order and Ordinary Climatological Stations 


In the previous sections, considerable mention was made in 
regards to the types of observing programs at different stations. This 
was due to the variations occurring in the characteristic curves among 


stations within reasonable proximity of each other. These variations 
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67 
or differences were most inherent among stations along the southern sec- 
tion of Alberta. In some instances, the discrepancies were explained 
through meteorology, in other cases, topography and: elevation were 
assumed to be factors, but in some situations it was presumed that 
errors in measurement or the type of observing was the key factor. 

With the definitions of the different types of stations given in 
Chapter IV, a comparative study can now be made between a first-order 
and a nearby ordinary climatological station. Consider the two local- 
ities of Medicine Hat and Seven Persons, separated by a distance of 15 
miles and having elevations of 2365 and 2480 feet above mean sea level, 
respectively. As indicated by Figure 32, the profile between the two 
localities appears smooth, thus topography should not be a major factor. 
It could be assumed that the distribution of precipitation should be 
the same for both stations. 

A comparison of the stations shown in Figure 48 indicate that 
a large variation exists between the two Ay curves. While both stations 
show the same seasonal variation, the Ay curve for Medicine Hat has a 
factor range of 1.2 to 2.0 greater than the hy curve for Seven Persons 
for the entire year. Due to the nature of the type of precipitation 
over Alberta, especially in the case of convective storms, one could 
expect some differences between nearby stations, say in the order of 10 
to 15 percent. With a variation as high as 50 percent, one can only 
assume that the nature of observing is the over-all factor in account- 
ing for such differences between Medicine Hat and Seven Persons. 
First-order stations are staffed with trained personnel providing contin- 


uous coverage of the meteorological elements while most ordinary clima- 
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68 
tological stations are maintained by voluntary observers who are not 
compelled to provide continuous observation. It therefore can be 
assumed that many ate providing small amounts of precipitation are 
not recorded at ordinary climatological stations because of wind-action, 
evaporation and other effects. This could be a partial explanation for 
the discrepancies between the nearby stations. 

An inspection of Figure 49 reveals that the Ay curve for Seven 
Persons has the same range of difference that was discovered in the 
comparison of the Ay curves. Referring to Table II, the mean annual 
precipitation for Medicine Hat and Seven Persons is 14.29 inches and 
13.79 inches respectively, but the total number of days with measure- 
able precipitation varies from 93 days at Medicine Hat to 66 days at 
Seven Persons. Rather than accept that more precipitation is expected 
per storm at Seven Persons because of the lower ro values, the author 
is compelled to believe that the observing programs at some ordinary 
climatological stations fail to report the total precipitation at their 
‘station 

Therefore, for ehe characteristic curves of an ordinary climatolo- 
gical station to be representative of the precipitation in its surround- 
ing area, and adjustment factor may be required. The magnitude of such 
a factor and where it can be applied would require extensive studies 
with an intricate precipitation network. But for practical purposes, 
such curves could be used for rough estimates. 

In the event that the comparative study of the characteristic 
curves for Medicine Hat and Seven Persons was an extreme situation and 


to avoid using Seven Persons as an example for all ordinary climatolo- 
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69 
gical stations, an inspection was made of the hy and do curves for Fort 
MacLeod and Lethbridge. The elevation chart of Figure 32 shows that 
there is a minor difference in elevation between the two stations and 
the topography is of a slight rolling nature. While there will be some 
effect of topography on the distribution of precipitation at Fort 
MacLeod because of its nearness to the foothills, one can assume that 
the distribution of precipitation at the two localities should have a 
high degree of similarity. 

Referring to Table II again, the mean annual precipitation for 
Lethbridge is 17.23 inches and for Fort MacLeod is 17.67 inches while 
the number of days with precipitation for the former is 110 days and 
only 88 days for the latter station. This difference in the number of 
stormy days is in agreement with Figure 50 which illustrates that the 
AL values for Fort MacLeod are less than in the case of Lethbridge 
except for the latter part of August and all of September. The great- 
est differences between the curves are during the winter months and it 
can be presumed that the high frequency of chinook winds in this area 
may have a considerable impact on snow measurement. The difference in 
values in the period of late July to early August can be attributed to 
the number of convective storms that are recorded, especially those 
storms producing small amounts of rain but are lost due to evaporation 
if times between observations are not less than 24 hours. It can also 
be seen from Figure 50 that the same seasonal variation occurs at both 
stations. 

Since Lethbridge and Fort MacLeod have approximately the same 


mean annual precipitation, one would expect that with Fort MacLeod 
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having the lower Ay values, it would have lower dy values than Lethbridge 
and this is generally confirmed by Figure 51 which gives the comparison 
of the do curves for both stations. Similar to the Ay curves of 
Figure 50, the greatest difference of the do values occur during the 
winter months when a greater part of the precipitation is snowfall. 

As can be seen from both aT and hy curves for Lethbridge and 
Fort MacLeod, making adjustments to the characteristic curves for an 
ordinary climatological station would be difficult without further study 
because of the time and space variability of precipitation. An im- 
portant factor that has arisen from the study of the curves for both 
Lethbridge and Fort MacLeod is that the measurement of snowfall may not 
be precise at ordinary climatological stations in or near mountainous 
regions. As a result, the amount of snowfall may be under-estimated due 
to the problem of measurement and this may have a serious effect on the 
conclusions drawn by the hydrologist or hydraulic engineer who makes 
full use of precipitation amounts in his calculations of river runoff. 

From the analysis of both sets of stations, it has been shown 
that the measurement of precipitation at ordinary climatological 
stations is inferior to that of first-order stations and a certain 
degree of caution should be exercised when employing precipitation 


records for important analytical considerations in hydrology. 


5.4. Effects of Elevation 

Besides its relationships with time and area, precipitation dis- 
plays variations with other parameters. These may be physiographic, 
like height, slope, aspect, and distance from the sea, or dependent on 


airmass properties. On occasions, rainfall has been related to one 











oe oF oe . _ 


~ 


4 . i de : 7 
. are : - : y 
| eehbedsged nada genipav .\ tewel oval biaeW NH ,eodtev — 


ie | 


Td -asvmuo 2A att te Jotindia dtod ~menian edt Bo! 
4 ee _ 


erig Boao ied b yuvDo0 Bu v id te pineisdTrbh Feeds ens. oe orgt 5 
t ~t eed 
¢ - : ¢ al in 


baiwonme ef mal MI O tre te tes15 ther 2a Dereey twinkw 
, - 


j E 
at { ! 78) ; 2 , . ~ 
ni ‘ a Ve 
r 5 . ee g t ‘ | 
- 4 4 t a - 
’ { ' » alt - 
’ 
< 7 b 4 Lic om 30 
% 
-_" fo oe nt 4 ~ @ LY 
e 5 i? & 2I6 ‘ 4 4 TOVSeZ 
j e 
o ” ‘ - 
” ’ + = *% 
® C f O ies ine ( 5 2hE Lb a0 Beit 
~ ° \ 
A r . , t « ”” ? 
‘ 
‘ « te a r we 
' . deat | 
ty, 
7 é s 
” 4 . 
l = = os 
{ 
“F - 1 a 
, ~ yi : ~ a 
" ry * 
» 4 a + ) 
c m —~>¢ + 
‘ { £3 ra ” f J 






enol tad 


seveiS } 


- de 


: Libs 
ptsaten 






\ 


a bus ats sly te 4 
7 ui 





eal 
particular variable but with limited success. As in most situations, 
there is more than one significant variable. For example, an empirical 
relationship can be established between rainfall and elevation, but its 
use is limited to a specific locality and cannot be generalized with- 
out involving some other physiographic variable. It has been shown by 
authors in meteorology (Wiesner, 1970) that there is an increase of 
precipitation with elevation. The main feature about rainfall, as a 
mountain range is ascended, is not the higher intensities received but 
the longer periods of rainfall, which are often of light intensity. 

The need to study the variation of precipitation with elevation 
in Alberta stems from the fact that most important watersheds are em- 
bedded in the mountainous regions. An important future concept would 
be to find the degree of variation of precipitation with height and the 
critical elevation where the maximum amount of rain or snow is likely 
to occur. Then this information could be related to the elevations of 
any watershed under study. Eventually, a representative distribution 
ef precipitation could be obtained without involving an extensive and 
expensive precipitation gauge network. 

In mountainous regions, the variations in precipitation make 
assessment particularly difficult especially because sites for rain- 
gauges are few because of inaccessibility and are often non-representa- 
tive. For this particular study, the problem of evaluating the effects 
of elevation on precipitation for the Alberta mountains was quite diffi- 
cult as eitewadthor wished to compare two aeat tous within the mountains 
that would have the following criteria: 


1. The stations be separated by a distance not greater than 10 
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72 
or 15 miles. 
‘2. The elevation difference being 1000 feet or greater. 

3. That at least 25 to 30 years of continuous precipitation 
records were available. 

4, The time period of the records be the same for both local- 
ities. 

5. That both stations were of the first-order type providing 
continuous observations. 

To find a series of localities in Alberta that would satisfy all 
the above criteria was impossible. The two stations that could satisfy 
at least three points of the criteria were those of Cardston and Carway. 
Some information about the stations is given in Tables I and II. The 
elevation contour map of the Cardston-Carway region illustrated by Fig- 
ure 52 shows that Cardston is exposed in a clockwise direction from 
northwest to southeast and that Carway is embedded in a set of hills 
with elevations between 4400 and 4500 feet above mean sea level. Ideal- 
ly, for a better precipitation-elevation study, the station of Carway 
should have been located on the ridge with elevations of 4750 feet and 
greater. The elevation chart of Figure 53 is a north-northeast by 
south-southwest line between the two townsites and it indicates that 
Carway is slightly protected by the ridge to those winds from the north 
and northeasterly directions. In the winter months, those wind direc- 
tions have a tendency to cause upslope orographic precipitation, espe- 
cially paves ue Thus more snow is expected to accumulate on that 
ridge as compared to Carway. 


In order to determine the effect of elevation on precipitation 
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between the chosen stations, the characteristic Ay and hy curves were 
produced and compared as given by Figures 54 and 55, respectively. It 
can be seen from Figure 54 that both stations experience the same season- 
al variation and, except for the time interval of February to May in- 
clusive, the Ay values are nearly identical. The lower values of aT 
for Carway for that time interval may be due to some meteorological 
phenomenon, such as prevaiitne flow of moist or dry air. By examining 
Figure 55 and making reference to Figure 54 and Table II, the following 
deductions can be made: — 

1. During the period of June to September inclusive, both 
Stations receive the same amount of pee at hier a 

2. For the month of May only, Cardston is expected to receive 
slightly more precipitation as it has a Wiener Ay value but the same 
Ao value as compared to Carway. 

3. For the remaining months of October to April, the differences 
in the ho values in telation to the differences in the aT values suggest 
that Carway receives more precipitation per storm which during that time 
interval would be mainly in the form of snow. 

The above information is confirmed by the mean annual values 
given in Table II which notes that at the townsite of Carway, the mean 
annual values of total precipitation and snowfall are higher, but that 
the mean annual total number of days with measureable precipitation is 
less than i the case of Cardston. 

White the elevation difference between the two localities is 
only 685 feet, the change in mean annual precipitation amounted to 


slightly more than two inches indicating that greater amounts of pre- 
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74 
cipitation can be expected with height or near mountain ranges. While 
this concept has been proven by many scientists in meteorology, the 
purpose of this section is to indicate how the time distribution of 
precipitation is affected by changes in elevation. This was accomplish- 
ed through the development of the characteristic curves which demon- 
strated, at least for this particular area, when the differences in 
precipitation amounts could be expected. 

It should be remembered that the concepts developed in this 
section may not be applied in general because of the space variations 
of precipitation and also that the two stations studied were of the 
ordinary climatological type. Because of the nature of observing at 
these stations, randon and systematic errors are possible and this may 
alter the situation to some degree. Nevertheless, those hydrologists 
wishing to develop a complete precipitation distribution over a 
mountainous watershed must consider the elevation factor and may have 


to be satisfied with any type of station. 


eds Number of Years of Record Required 


Realizing that this study would involve the statistical treat- 
ment of daily precipitation for a number of stations, it was felt that 
at least 30 years of records would be required as a larger sample size 
would provide a smoother representation of the variable in question. 
For those stations that contain records of considerable length, a reason- 
ably smooth set of characteristic parameters could be obtained. But 
bastres te many stations (including Alberta) where the length of record 


is short and may not be continuous during some periods. For those sta- 





» 
4 = 7 
iy 4 r aa 4 pon ° i! : sn¢ of b>. L, 7 45 was ad 
~ : ; - 
f, i 
13 Bi Di racy A yory oad 
\e ; / 
f r ! ha ‘| 7 ’ ba 
4 . i om We a 
a2 2 } f f 
: ye 
ed 
4 5 adh 4 
' 
, 
‘ 
i 
‘ ef 
» 4 _ 
‘ 4 ‘ > hig 
? 
, 
~ i Pr ftis iS 
] : 
, 
' 
( 
: ‘ bead * 
i 
{ w 
. ' 4 Ls | “ 
~ if C { Ji¢ y| 
i . 
~— - ¢ —— 4 
34 : = - uy ro 
é 
7 ' te 5 p 
AGsIesGp Ls F Lieinsseties +91 . 
: ; b 
ROBBY 2 .doenel: 4Tdeysbienos 26 abeases ft 
RONSS1 BH ~AIBIOL slosyablano> 7 PTOIG4 rke wrod sci anotiesa ‘pac 
7 yy 7 4 — war 





pice 
> * a “a1 7 
bate Bydd rey Bleos sxate neq Signet tedon tans to, di me — w 
wear z 
OME ES iss: 
to: tages! in sian t tredih gaxtetomk) | 















Seons sou iaq & 
_ a ras 
i ar —_ Tu . 

§ 





5 
tions, some suitable length of record would have to es determined which 
could cause a variable to be representative or characteristic of that 
station. 

To provide some insight into the question of the length of 
record required to obtain a suitable graph of Ay or do versus time of 
the year, the record for Edmonton was divided into two sections; one of 
ten years from 1959 to 1968 and one of fifteen years from 1954 to 1968. 
The values of Ay were computed for both periods and compared against 
the 3l-year time average while the ro values were computed for the fif- 
teen-year period and also checked against the longer period of record. 
The comparisons are plotted on Figures 56, 57 and 58. 

An examination of Figures 56 and 57 shows a close similarity 
between the ten-year and fifteen-year time average curves with the 3l- 
year average curve. However, the fifteen-year curve has a closer and 
smoother fit to the longer period curve than that of the ten-year record 
length. Referring to Figure 13, it is apparent that the shorter periods 
of record show a greater scatter about the time average curve. Several 
important features are evident with the inspection of Figures 13, 56 
and 57 and some are as follows: 

1. The similarity between the Ay graphs for the periods 1938- 
1968, 1954-1968 and 1959-1968 indicate that no significant change in 
the number of measureable days with precipitation has occurred since 
1938 and it may be safe to assume that little change has occurred since 
observations began in 1880. | 

Bie nieeid no doubt that a very long period of record would 


eventually provide a reasonably smooth curve. By inspecting Figure 57 
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76 
and noting the close similarity between the curves, it is quite prac- 
tical to say that the averaging of data for a shorter period of record 
will provide a curve that is a close approximation of the smooth curve 
that would have been obtained if an infinite period of record were 
available. 

3. With the concept of using shorter periods of record for the 
development of characteristic curves for different variables associated 
with the precipitation phenomenon, one can now perform analysis on 
newer stations that may be located in remote areas and therefore allow- 
ing a greater coverage for any specific zone or region. 

In Figure 58, the characteristic do curves for Edmonton are 
shown as the smoothed curve of the computed points for the period 1954- 
1968 and the time average curve that was drawn in Figure 26. Both 
curves have the same seasonal distribution and can be classed as being 
similar. For the period 1954-1968, the do values except for the months 
of July to October inclusive are higher indicating a smaller amount of 
precipitation per storm for the most recent period. Since the precipi- 
tation gauge for Edmonton has not been moved any considerable distance 
in the past fifteen years and since there is a very good agreement of 
the Ay curves from Figure 57, it can be assumed that no major change in 
the climate has occurred. However, the possibility of a small fluctua- 
tion or oscillation in the climate should not be overlooked. 

It is unlikely that the difference in the amount of precipitation 
is purely Shee of chance and some other factors must be the cause 
for this change. It would be difficult to speculate as to what these 


factors are, but one could take into account the rapid expansion that 
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7 
the City of Edmonton has experienced in the past fifteen years or the 
increased number of pollutants that exist in the airmass as of recent 
times. As was noted with the comparison of the Ay curves for Edmonton 
for the different time periods, there is a greater scatter about the 
mean value curve of the Ay values for shorter periods of record. In 
order to establish the effect of the differences in dy between the 
periods of 1954-1968 and 1938-1968 on the expected amount of precipi- 
tation [E (x, )] for Edmonton for any time interval, it was decided to 
divide the records into two parts and compare the resulting curves 
against those of the longest period. 

Consequently, the records were divided into the time periods of 
1938-1953 and 1954-1968 and compared with the total length of record 
consisting of the years 1938-1968. The characteristic Ay and Ay curves 
were drawn on Figures 57 and 58. However, the plotted points apply 
only to the interval 1954-1968. It can be seen from Figures 57 and 58 
that while all curves generally have the same seasonal variation, there 
are differences and similarities during different portions of the year 
for all three time periods. Most noticeable is the departure in the do 
values of Figure 58 during the first six months of the year and then a 
good similarity in the values from July to October and finally some 
small differences in November and December. It would appear that the 
smaller amounts of precipitation per storm during 1954-1968 as compared 
to the previous sixteen years of record for the first half of the year 
would suggest that possibly a small oscillation in the climate had 
occurred over Edmonton in the past 30 years. A verification of such an 


oscillation would require further study of other meteorological param- 
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78 
eters. 
To determine the degree of error that would result in the expect- 
ed amount of precipitation for Edmonton for any time interval using a 
shorter length of record, the months of June and September were selected 


and the data presented in Table VI. 


TABLE VI 


EXPECTED AMOUNT OF PRECIPITATION FOR EDMONTON FOR JUNE AND SEPTEMBER 


Per iod 
Month ui 2 3 
1938-1953 1954-1968 1938-1968 
June 2.03 2.47 2.87%* 
September 0.96 ree Ls43 


** Amounts are expressed in inches. 


The amounts shown in Table VI were obtained from the Ay and ho values 
extracted from the curves for the three time periods. The percentage 
error of periods 1 and 2 with period 3 for the month of June is 5.6 
and 13.7 percent, respectively. For the month of September. the error 
of periods 1 and 2 with period 3 is 15.1 and 9.5 percent, respectively. 
The different values demonstrate that significant deviations occur in 
the expected amount of precipitation when shorter periods of record are 
used. If necessary, certain adjustment factors may have to be devel- 


oped and applied to the analyzed data of records of shorter length. 


5.6 Summary 


With the development of the characteristic I and A, curves for 
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79 
the precipitation variable for any station, some knowledge can be 
attained as to how the precipitation is distributed in time. By com- 
paring stations within a certain zone or region, the d, and do curves 
can give some information about the distribution of precipitation in 
space. Grouping of stations that contain some similarity and show the 
same seasonal variations allows the classification of geographic zones 
and therefore a valuable concept in the time and space distribution of 
precipitation over an extensive area. Further study by the hydrologist 
of the precipitation phenomenon with other meteorological parameters 
will enable him to determine the overall effect of the distribution of 
precipitation on the course of any river from its watershed to the 
eventual termination in any large body of water. The characteristic 
curves can also be useful in determining, to some degree, the type of 
precipitation that occurs for any specific locality. 

Any project study involving the accurate determination of run- 
off using precipitation data will require that these data be reasonably 
accurate. Such accuracy can only be attained if the station reporting 
precipitation has quality control, that is, the station is staffed by 
trained personnel who provide a continuous observation using the most 
recent and sophisticated instrumentation. Of the several types of sta- 
tions in Alberta, the most common are the first-order and ordinary 
climatological stations. Comparative studies of the above stations 
when in close proximity of each other reveals that some departures do 
exist and ne Cases were quite excessive. This is mainly due to 
the differences in observational techniques, therefore, some caution 


should be heeded when data are being used from ordinary climatological 
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stations. It has been shown that some of these stations are likely to 
report an undercatch of the true precipitation. Whenever possible, 
some correction factor should be applied to ordinary climatological 
Stations to make their records compatible with first-order stations. 

Knowledge of precipitation in mountainous terrain, especially 
over major watersheds, is vital in the determination of e595 runoff 
for flood forecast ing and other purposes. It is generally accepted 
that precipitation amounts per storm increase with height to a critical 
elevation and then decrease above that level. The study of the effects 
of elevation on precipitation becomes important in those remote areas 
where the hydrologist has only a few stations available, but by 
account ing for the changes in elevation, a reasonable distribution of 
precipitation can be obtained. 

The determination of che minimum number of years of record 
required ba give reasonably smooth Ay and do curves was undertaken as 
many stations in Alberta have short lengths of continuous records. It 
has been shown that a record length of fifteen years is satisfactory 
in providing characteristic curves that would be representative of a 
particular station. Although appreciable deviations did occur when a 
shorter period of record was used as compared to a longer length of 
data, an adjustment factor may be introduced or an improved method of 
smoothing could be initiated. It should be remembered that all of the 
smoothing of the curves in this project was through the use of moving 
means and then a final correction by eyesight. 

And finally, it should be noted that while some meteorological 


theory was used to explain anomalies between stations, especially in 
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mountainous regions, these same anomalies may have been the result of 
the type of observing program carried out at the station in question. 
It is impossible at this time to say what factor contributed most to 
the anomalies, but it is hoped that further research into the meteoro- 
logical aspects of precipitation and as well the amount of undercatch 
of precipitation at ordinary climatological stations will provide more 
concrete evidence as to the behavior of precipitation distributions in 


Alberta. 
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CHAPTER VI 
COMPARISON BETWEEN OBSERVED AND THEORETICAL DISTRIBUTIONS 


6.1 Introduction 

Further to the development of the characteristic Ay and Ay curves 
based on the techniques illustrated in Chapter IV and with the extended 
discussion of the curves in Chapter V, it will be the purpose of this 
chapter to determine the theoretical distributions for certain eta one 
by the substitution of Ay and Ao values in the theoretical equations 
derived in Chapter II. The theoretically derived distributions will 
be compared to the actual distributions and then it can be determined 
whether the hy and ho values extracted from the average curves will 
provide distributions that are in close proximity to the population 
distribution compared to the case where the actual computed values of 


A, and Ao were used. In some instances, especially over short inter- 


iy 
vals of time, the scatter about the average curve may be appreciable. 
Thus one should expect noticeable discrepancies between observed and 
theoretical distributions. In the tables referred to in this chapter, 
for simplicity, abbreviations will be used; 'obsvd' representing the 
observed data and 'theor' denoting the theoretical distributions based 
on Ay and ro values derived from the characteristic curves. 


In the succeeding sections, comparisons will be made for single 


stations, groups of stations, for differences between nearby stations 
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and it will be investigated whether or not the theoretical distributions 
can be used to determine the precipitation regime over geographic areas. 
The ultimate goal of such a study is that of providing a first estimate 
of precipitation behaviour over a region that has little or no precipi- 
tation and hydrological instrumentation. This may be accomplished 
through a certain method of transposition of characteristic curves from 
a nearby station and accommodated with the proper adjustment factors 
which might include elevation, topography, orientation, hemispheric 
location and some knowledge of the climate. 

The tables providing the comparisons of the observed and 
theoretical distributions and which are referred to in this chapter 


are given in Appendix I. 


622 Number of Storms During a Specified Time Period 


From the examination of the graphs illustrated by Figures 6 to 
18, it was discovered that the function A, (t) computed from the distri- 
Bre ton of E(n,) was not constant with respect to time. With reference 
to the general equation for the mathematical expectation of the random 


variable Ne (Equation 2.2), 
EM) = f'n, (ddr 
Soe Gee al 


it appears that, since E(n,) is merely the area under the characteristic 
Ay curve over the time interval considered, then the time average value 
Q,) of Ay over the time interval can be employed. In other words, 
the characteristic curves represent the time average values of hie 


Therefore, 
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NM )drem A (Heer e) 
a A ey ib Pas a 


and which allows the use of the Poisson distribution with parameter 
Xt GE to is taken as zero) for the distribution of nee 

To evaluate the effectiveness of the Poisson distribution in 
describing the precipitation phenomena, the observed frequencies and 
the ‘ARS frequencies were computed for the first five days of 
each month for Edmonton as given by Table VII and for the first five 
days of January for the four stations along the central Alberta zone 
as shown by Table VIII. From Figures 6 to 18, it is noted Ae a 
seasonal variation persists at every station and in most cases a sudden 
change occurs in the Ay values between May and July. In order to test 
the applicability of the theoretical equation during the time interval 
of a sudden change in Ay observed and theoretical frequencies were 
computed for a 20-day period. For Edmonton and Edson, this period was 
from June 1 to 20, while May 21 to June 10 was selected for the stations 
of Entrance and Jasper. The results are shown in Tables IX and X, 
respectively. In all cases, the Ay values (shown in the tables) are 
average values over the time interval from the characteristic Ay curves 
of Figures 13, 14, 15 and 16. 

An overall inspection of Table VII reveals that there is a 
reasonable agreement between the observed and theoretical distributions. 
The areas of disagreement are those for the frequency of occurrence for 
one day where the theoretical equation often gives an over-estimation 
of occurrence and in the case for two days, the theoretical distribution 


shows an under-estimation of the number of occurrences. It is presumed 












a ¥ ae ae 
_ 
a 
4 
ed = | 
1 ~ ? 
~~ : 
, i) p ~ 4. oe = 
(3 ~—o) 4 © Shir, \ ‘ | 
; oO I o a? 7 
4 
5 ; = 
} 7 i 7 : nig ahs | ar dv, u wit Tees 
' 
: ' 0 fe »2 ( sc 
Me 
1 e ; eee! 4 5 
j 4 ‘ Jan aya Mf 
‘ i 
: : . 
4 ‘ 
> +8 a fe jac 0) 
i j a3 i r 
oe 5 ' 4 — : an 
re j > a7? “ rp A & LY 4 iqgn sta 
: a 


mgMn « ab “9y : nahfos ® 36 


er 7 . ; mils “ine batt : «! i » wot to suqao> 
_ 
pst! , ) ‘wl: © aM 972 f o1 J boul work ; 

| ihe 

Lua BAY taqasL ae caz3al ‘to 

* ae - 

(a ) + i ( od? 4a ile st ew wine - 













7 a . . : ; : " ~% 
penolaudtsjarb Lmorisionis ina bevageds ais naw? @ James fgs 
> ; ’ 7 *'s ae 
— - ee ol sofeyws sa 70 bared vee) ‘ i SA3 07 sidaaetadl 376 ses Fak a0 
a wr, = : we 
er ~y Pe - 


a ; oe 
AB potaynasio « Lg Beer od 39d ae 


us wT 
est ai mt a. ws ih a — © ale ih 23 yee ho 





85 

that, by choosing the actual Ay value computed from the distribution 
of E(n,) instead of the Ay value selected from the characteristic 
curve, there will not be much change in the differences between the 
observed and theoretical distributions. 

The reader can recall from previous chapters that when daily 
precipitation data was used, there is the inevitable over-estimation 
of the number of storms when weather systems persist for more than one 
day. It is often the case in the analytical treatment of precipitation 
data that when two or more continuous days with measureable precipi- 
tation were counted, it was suspected that one major storm was occur- 
ring at the time. This type of situation is very prevalent over western 
Canada where well-developed weather systems move slowly eastward during 
certain seasons of the year and provide several days of precipitation 
to any specific locality. This condition has some bearing on the 
assumption made in Chapter II which stated that the probability a storm 
will occur in the interval t+At does not depend on the number of storms 
up to time t. This assumption of independence was required so that the 
theoretical equations could be developed (Verschuren, 1968; Todorovic, 
1967). With any weather system lasting for several days, it can be 
accepted that the probability of rain in the interval t+At is dependent 
on whether or not it rained up to time t. In other words then, the 
Poisson distribution equations developed in Chapter II based on the 
assumptions must be considered as an approximation only if daily pre- 
cipitation etext are used. As a result, differences in the observed 
and theoretical distributions should be expected. 


Any improvement of the theoretical distribution equations to 
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provide results similar to the actual or observed distributions will 
require the incorporation of meteorological theory involving precipi- 
tation regimes. Referring to Table VII again and accepting that the 
departures between observed and theoretical frequencies are small, 
then one can consider at least for the station of Edmonton that the 
theoretical distribution is a rather good approximation of the number 
of days with precipitation during the first five-day period in each 
month. It is certain that the theoretical Be oas would provide 
reasonable results for the remainder of the year. 

To determine whether the theoretical distribution equations are 
a good approximation for the other stations, the Ay values for the 
interval January 1-5 were extracted from the characteristic curves for 
‘pasedt Entrance and Jasper, and applied to the proper equations. The 
month of January was chosen because the oF value for all four stations 
displayed a constant value as compared to the other months. The 
results are given in Table VIII. It is noted that for the higher A, 
values, for instance Edmonton and Jasper, there is a good agreement 
between observed and theoretical frequencies; but for the lower values 
of Ay for Edson and Entrance, there is a poor agreement, especially 
for days of zero precipitation and for one-day occurrences. The 
discrepancies can be accounted for by the fact that since precipitation 
gauges often to Fail to record small amounts of precipitation, it must 
be expected that, when observed and theoretical distributions are com- 
pared, enewoneeheeal distribution will generally show a larger number 
of days with no precipitation. 


Mentioned previously was the problem of'a varying Ay dur ing 
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marked changes in the séasonal variations and what effect it would have 
on the theoretical distribution. Table IX and X give the results for 
the 20-day interval that included the maximum change in ry for the four 
central Alberta stations. The oF value used is the average value of 
the characteristic curve of the selected time intervals and should be 
a good approximation of the actual AG value based on the distribution 
of E(n)> therefore a reasonable agreement should be expected between 
the observed and theoretical frequencies. Inspection of Tables IX and 
X indicates that while differences in frequencies do occur, the depar- 
tures are not serious and it can be accepted that the use of an average 
value of aT will be representative for those Bie, Bee anaes there is 
a marked seasonal variation. It is possible that a longer period of 
record would show a better correspondence between observed and theoret- 
ical distributions for a Ay curve that has marked seasonal variations. 

Another interesting test is that of evaluating the effectiveness 
of the theoretical distribution for longer intervals of time. Thus 
Tables XI and XII were developed to compare the observed and theoret- 
ical distributions for the central Alberta stations during the 30-day 
period of January 1-30. The Ay values shown in the tables were ex- 
tracted from the characteristic curves and should be similar to the 
averaged Ay values obtained from the observed data. This suggests that 
discrepancies between observed and theoretical frequencies should be 
minimized. 

In the situations for Edmonton and Jasper, both having high ArT 
values, there is a reasonable agreement between the observed and 


theoretical frequencies. In the case of Edson and Entrance where lower 
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88 
Ay values were experienced, there is disagreement between the two 
distributions. The differences in the latter case might be attributed 
to the following factors: 

1. The Poisson distribution equations are ineffective at low hy 
values. 

2. The observing program at Edson (a first-order station) under- 
estimated the number of days with measureable precipitation, especially 
snowfall during January. Actually, Edson did not become a first-order 
station with full instrumentation until February of 1960 and therefore 
a good rosaibility exists that during the period of record 1929-1959, 
many stormy days with small prec ipitation amounts could have been 
neglected. 

3. With Entrance being an ordinary climatological station and 
since it has been shown that such stations are likely to under-estimate 
the number of stormy days, then a constant error over the years of 
record may have some significant effect on the observed frequency. 

4. The larger departures between the observed and theoretical 
distributions are for zero and one-day occurrences of precipitation. 
The theoretical equations based on their assumptions are unable to 
account for the anomalies. 

Discounting point 1 as unlikely and considering that points 2 
and 3 have some effect, the emphasis should ferpleced on point 4 to 
account for the differences in the two distributions. There is a good 
likelihood that the probability of rain or no rain on any day has its 
greatest impact on those days with zero precipitation and on one day 


occurrences. 
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89 
In summary then, it should be accepted that the Poisson equation 
provides a theoretical distribution that is a good approximation to the 
precipitation regime for an area with the condition that the character- 
istic curves developed for the area are representative of the precipi- 


tation variable. 


6.3 Amount of Precipitation for a Given Number of Storms 


The analytical development of ry and do values, the graphical 
representation of the characteristic curves, and the introduction of 
theoretical distributions to approximate precipitation behaviour in the 
previous portions of this study has now led to the most important con- 
cept of hydrological and water resource planning. It is the determina- 
tion of the probability of occurrence of a quantity of precipitation 
less than or equal to a fixed amount in a selected time interval during 
a particular time of the year. The important application of the theory 
presented in previous chapters will be in the situations where some 
knowledge will be required about the amount of precipitation to be 
expected in an area with no instrumentation or any precipitation 
records. This could be accomplished through the transposition of 
characteristic curves of a nearby area and these curves would provide 
the parameters necessary to determine the probability of an expected 
amount of precipitation. | 

The theoretical distributions of the random variable X. for the 
four central Alberta stations were computed with the aid of a computer 
program and the resultant distributions were compared to the observed 


values. The month of January was selected as the LT and ro values for 
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90 
all the stations were reasonably po detane in comparison to the remain- 
ing months. Thus t, was chosen as January lst., 0 hours. For this 
particular month, the average values of hy and Ao selected from the 


characteristic curves are as follows: 


2G AD) 
Edmonton 385 v2 
Ed son ay sity 6.35 
Entrance Ap Ahe) 8.00 
Jasper 2300 SU 


The results of the calculations and the comparisons are given in 
Tables XIII to XVIII where Xx, was considered for the time intervals 
t = 5, 10, 15, 20, 25 and 30 days. Basically, there is a fair agree- 
ment between the observed and theoretical eae for “hl the time 
intervals selected. A detailed study of the Tables XIII to XVIII 
reveals the following properties: 

1. For Edmonton, the agreement between the observed and theoret- 
ical frequencies is good for all time intervals which demonstrates the 
importance of quality observations from a first-order station that has 
proper instrumentation. 

2. For the remaining three stations, there is reasonable agree- 
ment during the time intervals of t = 5 and 10 days. For the longer 
time intervals & t (t = 15, 20 days, etc.), differences between 
observed and theoretical frequencies became more noticeable. Sugges- 
tions for the departures in the cases of Edson and Entrance were given 


in the previous section. An explanation for the departures for the 
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Station of Jasper is difficult, but it may be speculated that the do 
value chosen for January is not penresent adv ee icr the entire month. 
Comparison of the frequencies for the remaining months for Jasper 
would be required to determine the overall representation of the dy 
curve for this station and the mountainous area in the vicinity. 

One more measure of testing the expected value of precipitation 
over a specified location is that of plotting the theoretical and 
observed values of E(X,) and Var (X,). For Figure 59, the stations of 
Medicine Hat and Edson were selected for the month of January. It can 
be seen that the average values of X are in close correspondence with 
the theoretical values and there is also a good agreement in the var- 
jances for the first 15-20 days of the month, However, the variance of 
the observed values become much greater as the time interval increases 
suggesting a greater dispersion about the mean value than the theory 


would indicate. 


6.4 Comparison Between Nearby Stations 


In chapter V, the comparison of first-order and ordinary cli- 
matological stations was initiated because of the noticeable differences 
in the characteristic curves between Medicine Hat and Seven Persons. 

The two stations are separated by a distance of 15 miles and are of 
different types. To dersdee examine the precipitation distribution of 
this particular area in Alberta and observe the behaviour of the 
theoretical distributions, the month of August was selected for the 
two stations as it was the only month which did not have any missing 


records for either station. The comparisons of frequencies are given 
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92 
in Tables XIX to XXIV. 

The prevalent features of the tables are as noted: 

1. For all time intervals, the observed frequencies of days with 
zero precipitation for Seven Persons are greater than Medicine Hat. 

2. With the exception of t = 5 and 10 days, generally there is 
a poor to fair agreement between observed and theoretical frequencies 
for both stations. Some of the disagreement, especially for longer 
time intervals, may be due to the fact that average values were being 
used for a varying Ay value. 

3. The theoretical distribution under-estimates the days of no 
precipitation especially for the station of Seven Persons. Considering 
that Seven Persons is an ordinary climatological station and subject 
to certain types of error because of its observing program, it would 
almost be feasible to use the theoretical distribution of Medicine Hat 
as being more representative of the conditions at this locality. 

4. In the case of Medicine Hat, as t increases, the theoretical 
distribution has a tendency to show a higher frequency for greater 
amount 8 of precipitation than does the observed distribution. 

The above arguments apply to the month of August only and should 
not be construed to be representative of these stations for the remain- 
ing months of the year. Time and space limitations of this thesis did 


not allow a complete inspection of a station for the entire year. 


6.5 Applicability of the Theoretical Distribution to Geographic Areas 


The concept of combining stations where some lie within the moun- 


tains, others on the lee side and the remaining stretched towards the 
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93 
plains region may be questionable. However, considering that the mean 
annual precio catiog (from Table I) for the central Alberta stations 
are similar to each other and the fact that the characteristic curves 
were also reasonably similar, it would seem practical to have such a 
combination of values to be representative of any large area. 

The construction of Tables XXV to XXX were by the following 
procedures: 

1. The observed average is the average of the frequencies of 
all four stations used in the central Alberta section. 

2. The theoretical average denoted with index 1 was evaluated 
through the use of dy and Ay values extracted from Figures 39 and 42, 
respectively. 

3. The theoretical average denoted with index 2 was taken from 
Tables XIII to XVIII in which the theoretical frequencies were averaged 
among the four stations. 

Although there is a slight under-est imation by the theoretical 
distribution for days with zero precipitation, one can accept that a 
reasonable correspondence exists between the two distributions. Also, 
with a high correlation between the two theoretical averages, it is 
apparent that the combined characteristic curves are representative for 
the area between Edmonton and the mountains. Therefore, it can be 
suggested that if characteristic curves in a region indicate a similar 
type of precipitation distributions, then such curves may be combined 
thereby Beige a single representative curve which will be applicable 


to the entire region and possibly toa greater extent. 
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94 
6.6 Test of the Ay. and A» Values with Published Meteorological Data 

In earlier sections, the AF and ho parameters as obtained from 
the characteristic curves were used to derive the theoretical distri- 
butions. To determine the representativeness of the characteristic 
curves compared to published data, the Edmonton and Seven Persons rec- 
ords were chosen. Only monthly average values of Ay and Ao could be 
used to compare with published meteorological data. 

The results are shown in Figures 60 to 63. For Edmonton, the 
comparison of the average monthly values of the number of days with 
measureable precipitation using the Ay parameter shows close agreement 
with the monthly averages obtained from published data (Canada Meteoro- 
logical Branch, 1967). The average monthly values of the amount of 
precipitation using the ro values shows reasonable agreement with the 
monthly averages obtained from published data. The discrepancies for 
the months of June, July and December may be due to the smoothing 
techniques used to develop the characteristic curves and which may have 
some effect on the values of Ay and ry for those months. At Seven 
Persons, there are some differences in the number of days with precipi- 
tation but the average amount of precipitation shows a good similarity. 

The meteorological published averages are usually based ona 25- 
to 30-year record length within the period from 1930 to 1960. In some 
instances, the ees are based on record lengths as low as ten years. 
In the case of Seven Persons, the monthly averages BE some precipitation 
parameters yee based on a period of record of less than ten years. In 
this thesis, all characteristic curves were derived from record lengths 


of 25 to 31 years within the period 1924-1968. The differences in the 



















— 


4 : ; ee vis 
anit ts yh x isven4 aM pe bat ty seit ‘yi itw BL str te? , i aa’ r ada bh ry 


Md 
= meee es oA tee at aly oe Th kane okie sane 
1 = 


‘ i 7 ; r aor * 
i - -? 
| | are io take aanotdned sab ieee 
ast? bediatde e« eaetousrey ,4 bee aah? perottaad 7193 a 
ad a aT, : ' = 
) os 7 : -_ 
: | | rer 
f 4-ty } j 58 nen >i 1495 of Heel ease wavs al jo ttei sates 
A at & ~ “ sd - . 7 © . v - ts oy oa 7 
~~ his, f 
. : e 
F - aaz ra 
at ry semis #f3 ' Sqe Vt 3220965 fee? mz » los 1180 ws 
‘ : i” 
4 7 . t z ’ — P = , r 
ia 7 Oe =8 Yer see) uns odd ,a7eb Sepaticw © he 1eq=a07 3S 
a ~ 
. 
- 
ae mw 
i P freon 4 Pa ronn, A < weg — 
s¢ bivgos ,i fne .f to eoviav egeteva Yidsnom yIn0 .aneots eTaw eam 
> 5 2 
: i 7 F t fobw ST 
. ‘ = » — 
, 
© 
a a G @ aw ’ £ ,et 
! <* 
4 | 
i+ j : ; : pee pee a ey br’ 5 
: : 
oa ¢ 2 : Hi : 2 udev i AO8 1 } P 4 e149 
J : 
4 , ~ 
i ) i * .* - 4, , be 4s 
i pe * s ses 7 ’ 
cme! 644 20 . ei SeEte at... (86@L G 
art Tormey , 
it | wie Ge a 
4 5 Ar « 
- ; bs . 4 Lac0 
| ~ 
{= * 
om | *, § i re 7 ii Rai, af i 
‘ r , 
\ t 4 6 4 7 _ (sv hy ; 
uw o ; 
: > \ r 
: . 5 wh eet ; ws a id Si: 
ve 
mq) ny As ' 737i! > s7Ta 
¥ (imie Anog ews i s%G to -Javoees #28 iV 
” 1 
c : 
ao ‘ - j fy, ¢ a’, 
=} B:o feasd Vil us ia hertetiowg Laat gol 
a - te 
os i y 
= oad a A aA ff . wy Pe | - » » , i tw Y~ wee os 
: _ Siete Ha eX OY i ou we — be Mer al Wik at Srn3 ira is] Ww hd Ree ote wIO 14 aaa: 
~ a * * 





3 es wot x8 luiaiel Ons ‘= heaped ete 6 nateve ed} | 690 
7 : aot i - a1) 7 fo 
skgtnesg see 10 esgic pas glia siti ms sent a Foe 
; = -_ 

: =) : 7 - 7 
7 = ~ - ae Chit : 

ot “fins we a 3c 10 Bxdnee aececide ‘asin 










ee Tin 


95 
period of record and the number of years within the period explains 
some of the discrepancies between published monthly values and values 


obtained from the Ay and do curves. 
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CHAPTER VII 
APPLICATIONS 


7.1 Introduction 

For the hydrologist and hydraulic engineer who may have con- 
sidered meteorology as simply "it rained" or "it did not rain", the dis- 
cussion of the study of such concepts as the number of storms during a 
specific time interval, the amount of precipitation per storm, charac- 
teristic curves, observational programs, theoretical distributions and 
others has probably prompted him to consider where such concepts would 
be of practical value. Having observed that the use of the Ay and ho 
parameters derived from representative curves provided a reasonable 
correspondence between observed and theoretical distributions, the 
specialist can then incorporate the theoretical functions (with certain 
probability limits) into the analytical treatment of the precipitation 
variable. Most often the engineer will be interested in precipitation 
amounts received at the ground. Thus the distribution of the random 
variable Xx (the amount of precipitation during a specified time inter- 
val) would be of first priority, therefore discussion will be mainly 


confined to the application of this randon variable. 


dine Irrigation 


In this area of study, several types of professionals would be 
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97 
concerned about the expected amount of precipitation at any specified 
area. For instance, the agronomist's interest would involve the 
amounts of water required for different types of crops because of the 
seasonal effects of evapotranspiration whereas the engineer and hydrolo- 
gist would have to incorporate the seasonal fluctuations of precipi- 
tation into the overall design of an irrigation system. To demonstrate 
how the distribution function of the random variable X. can be used as 
a guide, consider a case of row-crop farming in the Medicine Hat region 
where the known minimum amount of water required to raise a crop will 
aid in determining the risk of crop failure without irrigation. 

Let it be assumed that the minimum amount of water required to 
raise corn in the first 20 days of August is equal to 2.00 inches. If 
the characteristic curves for Medicine Hat are representative for the 
area surrounding the city, then the aT and ro values are 0.246 and 
D240. respectively. With these values and t = 20 days substituted into 
Equation 2.10, the probability that x = 2.00 inches or less is .95. 
That is to say, 95 percent of the time, irrigation water must be sup- 
plied in order to maintain a crop of corn that requires a minimum of 
2.00 inches during the first 20 days of August. 

For the situation described above, let it be presumed that 1.60 
inches can be supplied by irrigation and that the remaining 0.40 inches 
must be available from precipitation during those 20 days. The distri- 
bution function for Xx, shows that crop failure ee iidenecnc 20 percent 
of the time if 0.40 inches is substituted for x in Equation 2.10. 
Therefore, from the above example, irrigation planning can be based on 


a certain loss-to-risk basis that would provide a cost-benefit ratio 
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suitable to all parties that are affected by an irrigation system. 


1e3 Precipitation - Runoff Relationship 


The importance of precipitation - runoff relationships is of 
concern to that engineer or hydrologist who is in charge of providing 
an accurate flood forecast especially in those regions where there are 
no hydraulic control structures and where river discharge data are inade- 
quate to determine peak flows in streams and rivers. River forecasts 
are also required for estimating inflow to reservoirs in order to 
permit the most efficient operation for flood control and other purposes. 
In addition, there is an increasing demand for day-to-day forecasts of 
river stages and discharges by those interested in navigation, water 
supply, stream pollution, and many other related fields. 

The tools of an engineer involved with forecasting include pre- 
cipitation - runoff relations, unit hydrographs, routing methods, re- 
cession curves and stage - discharge relations. Because of the im- 
portance of the time factor, great stress must be placed on the develop- 
Bene of forecast procedures that will enable flood warnings to be issued 
at the earliest possible time. Any forecasting service is dependent on 
adequate data. The development of the river - forecasting procedures 
requires historical data, while the preparation of operational fore- 
casts requires sufficient current information. 

Most often a design storm, commonly referred to as the probable 
maximum storm, is placed ayer a drainage basin and Brter losses due to 
infiltration, seepage and other causes are accounted for, the excess 


rainfall is converted to river discharge at the concentration point. 
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In other words, the probable maximum flood can be obtained from the 
probable maximum precipitation by subtracting loss rates and adding the 
snow melt with due regard to relationships with time. 

Judgements must be made in interpreting data and in reaching 
estimates of the probable maximum precipitation as currently there are 
no set of rules, graphs, and procedures whereby one can proceed step by 
step and necessarily derive an acceptable estimate. The application of 
the following methods may provide a reasonable possible estimate and 
they are: 

1. The storm model approach. 

2, The maximization and transposition of actual storms. 

3. The use of generalized data or of maximized depth, duration, 
area data from storms. These are derived from convective or orographic 
storms. 

4. The use of empirical formulae determined from maximum depth, 
duration, area data, or from theory. 

5. The use of empirical relationships between the variables in 
particular valleys. 

6. The statistical anlyses of extreme rainfalls. 

Unless the hydrologist or engineer is highly specialized in 
meteorology and has a capacity to apply some of the methods mentioned 
above, reliance must be placed on other statistical procedures that 
afford a shorter time in analytical treatment. It has been the purpose 
of this Pear to suggest that a reasonable degree of certainty can be 
obtained through the application of the characteristic Ay and Ao curves 


in combination with the distribution function for X when the statisti- 
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cal comparison of different design storms is necessary. 

The engineer who feels that depth, duration and area are the 
major factors in which a design storm would be used in determining a 
maximum probable flood, then the probability that the precipitation 
will be equal to or less than some amount x can be obtained through the 
use of Equation 2.10. The value of the depth of precipitation (x) and 
the duration (t) are the substitution parameters selected from any area- 
depth-duration curve. For example, if t is chosen as 5 days and the 
curve indicates a value for x of 0.50 inches, then the probability of 
this amount of precipitation being exceeded is equal to l - FL Gx), or 
Oc) Seif the AT and do values are 0.246 and 5.40, respectively. An 
important note to be remembered is that, since the values vary because 
of seasonal variations, the probability of occurrence of the selected 
design storm will not always be the same during different times of the 


year. 


7.4 Weather Modification 
The importance of the knowledge of precipitation distributions 
in the area of weather modification should never be neglected because 
of the future implications resulting from the modification of weather 
systems. Before any modification program is deemed successful or not, 
many facts must be searched by the scientists who are attempting to 
induce more precipitation from the cloud structure of a storm system. 
It can readily be seen that if the characteristic curves are 
accepted to be representative for an area or region for a considerable 


length of time, then the two variables: 
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1. the number of storms during a specified time period, and 
2. the amount of precipitation per storm, 
are going to assist the scientists in determining both short and long 
range effects of any precipitation modification program. The character- 
istic curves will have the capacity of directing pork the time and 
spatial requirements of precipitation modification, once it has been 


proven to be effective and feasible. 
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CHAPTER VIII 
CONCLUSIONS AND RECOMMENDATIONS 


8.1 Conclusions 

The sequence of this thesis, commencing with the development of 
the theoretical equations through the modification of the Poisson equa- 
tion, the derivation and discussion of characteristic parameters, and 
then the application of the theoretical distributions to describe the 
precipitation phenomenon, has provided the following conclusions: 

l. The assumptions used to derive the theoretical equations 
were not completely valid, but this problem is not serious and it is 
reasonabie to expect the theoretical distributions to be a good first 
approximation of the behaviour of precipitation. 

2. The two most important variables required to make the 
theoretical distributions applicable are: 

a. The number of storms during a specified time interval. 
b. The amount of precipitation per storm. 
The two variables are subject to seasonal fluctuations with a periodi- 
city of ane year and are represented by two characteristic curves that 
can be used to obtain an estimate of the probability of occurrence of 
an amount of precipitation in any time period. 
Se The attempt to describe the anomalies between characteristic 


curves for different stations with meteorology should not be accepted 
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CHAPTER VIII 


CONCLUSIONS AND RECOMMENDATIONS 


8.1 Conclusions 

The sequence of this thesis, commencing with the development of 
the theoretical equations through the modification of the Poisson equa- 
tion, the derivation and discussion of characteristic parameters, and 
then the application of the theoretical distributions to describe the 
precipitation phenomenon, has provided the following conclusions: 

1. The assumptions used to derive the theoretical equations 
were not completely valid, but this problem is not serious and it is 
reasonable to expect the theoretical distributions to be a good first 
approximation of the behaviour of precipitation. 

2. The two most important variables required to make the 
theoretical distributions applicable are; 

a. The number of storms during a specified time interval. 
b. The amount of precipitation per storm. 
The two variables are subject to seasonal fluctuations with a periodi- 
city of one year and are represented by two characteristic curves that 
can be used to obtain an estimate of the probability of occurrence of 
an amount of precipitation in any time period. 
3. The attempt to describe the anomalies between characteristic 


curves for different stations with meteorology should not be accepted 
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as conclusive as these anomalies can also be the result of the type of 
observing program at each station. There was no attempt to determine 
which factor had the greater influence. 

4. Although longer periods of record lend themselves to better 
statistical treatment, a minimum of 15 years of record in conjunction 
with proper smoothing techniques can be used to develop characteristic 
curves that would be representative of a longer period of record. 

5. The discrepancies resulting fn the comparison of the 
observed and theoretical distributions when daily precipitation data 
are used is greater for ordinary climatological stations than in the 
case of first-order stations. | 

6. The theoretical distribution, in most cases, under-estimated 
the number of days with zero precipitation and also over-estimated the 


larger precipitation amounts for a given number of storms. 


8.2 Recommendations for Further Research 

Due to the time limitations and the unavailability of more data, 
Mrerat problem areas of precipitation behaviour and regime could not 
be explored. Some of the more intersting ones are as follows: 

1. To improve the assumptions used in the derivation of the 
theoretical equations where daily precipitation is employed, the depend- 
ence of rain in the interval t+At on whether or not it rained up to 
time t should be considered. Also, a smaller time interval of precipi- 
tation measurement could be used. The use of hourly Been records 
would alleviate the problem in the classification of "storm" and 


"stormy days". Such a procedure would require the need of continuous 
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104 
recording devices as well as a properly chosen time interval that would 
define a storm or the break between storms. | 
2. In some applications of hydraulic engineering practices, 
there is a requirement for the ensues of the time period for a 
specific were: of precipitation to occur. Some examples are: 
a. Storage capacity for water management. 
b. Reservoir-filling time behind large hydraulic structures. 
c. Irrigation and agricultural crop practices. 

The determination of this variable requires continuous data. 

3. When the functions dj (t) and A, (x) were assumed constant, 
the theoretical distributions were based on Ay and Ay values that were 
linear with respect to time. This is not accurate when one considers 
that these functions are subject to seasonal variations. To accommo- 
date those situations where Ay and ro vary rapidly, the theoretical 
equations could be modified to follow a non-linear relationship with 
time and, therefore, provide a better agreement between observed and 
theoretical distributions. 

4, The involvement of meteorological theory of airmass circu- 
lation and storm movements would be an important asset in the stochastic 
treatment of the distribution of precipitation as it-would provide a 


means of establishing a classification system for precipitation. 
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APPENDIX A 


Program for Converting Raw Data for First-order Stations 
ee ee ee om ene OTe 
FORTRAN 


DIMENSION STRING(20) 
DATA RAIN/' — T'/ 
DATA ZERO/'0000'/ 

100 FORMAT(20A4) 

101 FORMAT(44H } ,4Al, 
132H ) | 
READ(4,100) STRING 
WRITE(7,100) STRING 

3 READ(4,101,END=2) PCPN 

IF(PCPN.EQ.RAIN) PCPN=ZERO 

WRITE(7,101) PCPN 

COMORS 

ENDFILE 7 

STOP 

END 


bo 
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APPENDIX B 


Program for Calculating Characteristic Curve Parameters for 


First-order Stations 


FORTRAN 


100 
101 
102 
#03 
104 
105 


gL 


No 


DIMENSION B(31,12,6),C(31,12,6),V(31) ,NAME(6) ,ND(12), Z(31) 
1, LMDA1(6) , LMDA2(6) 

REAL LMDA1, LMDA2 

DNTASND/O1m28451430/31,50, 31.31.30, 31,30, 31/ 

FORMAT('4',40X,'PROGRAM FOR FIRST ORDER STATION'//) 

FORMAT(17,6X,6A4) 

FORMAT (44X, F4.2) 

FORMAT (44X,17,4X,6A4,//) 

FORMAT (14,6 (2X, F6. 3, 2X, F6. 3)) 

FORMAT ('MNTH' ,2X, 'LMDA1',2X, 'LMDA2',4X, 'LMDA1', 2X, 'LMDA2', 
14X,"LMDA1',2X, 'LMDA2',4X, 'LMDA1, 2X, 'LMDA2' ,4X, 'LMDA1', 2X, 
2'IMDA2',4X, 'LMDA1', 2X, 'LMDA2', /) 

WRITE (6,100) 

READ(7,101) IP, (NAME(I),I=1,6) 

WRITE(6,103) IP, (NAME(I),I=1,6) 

WRITE (6,105) 

IFY=1938 

DO 1 I=1,31 

DO? Jal a2 

DO 11 M=1,31 

Z(M) =0 

KMAX=ND (J) 

IF(J.NE.2) GO TO 4 

IF(MOD(IFY+I-1,4).EQ.0) KMAX=KMAX+1 

DO 3 K=1,KMAX 

READ(7,102,END=8) V(K) 

IF(V(K).GT.9.99) GO TO 4 

IF(V(K).GT.0.0) Z(K)=1 

CONTINUE 

JMAX=KMAX/5 

DO 5S K=1,JMAX 

SUMB=0.0 

SUMC=0.0 

DO 6 L=1,5 

SUMB=SUMB+V (5 * (K-1) +L) 

SUMC=SUMC+Z(5* (K-1)+L) 

CONTINUE 

B(I,J,K)=SUMB 

C(I,J,K)=SUMC 

CONTINUE 

CONTINUE 

CONTINUE 
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APPENDIX B (CONT'D) 


DO 7 J=1,12 
IF(J.NE.2) JMAX=6 
IF(J.EQ.2) JMAX=5 

DO 9 K=1,JMAX 
SUML1=0.0 

SUML2=0.0 

DO 10 I=1,31 
SUML1=SUML1+B(I,J,K) 
SUML2=SUML2+C(U, J, K) 
CONTINUE 

LMDA1 (K)=SUML2/ (5*31) 
LMDA2 (K) =SUML2/SUML1 
CONTINUE 

WRITE(6,104) J, ((LMDA1(K) , LMDA2(K)) ,K=1,JMAX) 
CONTINUE 

STOP 

END 
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APPENDIX C 


Computer Output for First-order Station 
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APPENDIX D 


Program for Converting Raw Data for Ordinary Climate Stations 
FORTRAN 


DIMENSION IY(6),1IM(6) ,PCPN(6) ,STRING(20) ,ND(12) 
Tipon11 (6-9) ,12 (6,8) 
DAWARND e3ilee2 S551). 50, 31, 30,51, 51), 50,31, 50, 31/ 
DADAMREGSIMBLANK/' - D1) t 9 Cho t  opet '/ 
DATA ZERO/'0000'/ 
DATA DNEG/'- 99'/ 
100  FORMAT(20A4) 
102 FORMAT(I7,12,12,8A4,A1,A4, 8A4) 
READ(4,100) STRING 
WRITE(7,100) STRING 
28 DO 10 K=1,5 
READ (44102 END=8) ID(K) , EY(K), IMCK) , (T1(K,1),1=1,9), 
JROPN(K)i(T2 (K, 1), I=1,8) 
JK=1 
LM=IM(K) 
LY=IY(K) 
N=ND(LM) 
IF(PCPN(K) .GT.9.99) GO TO 28 
IF(LM.NE.2) GO TO 10 
IF (MOD((LY+1900) ,4).EQ.0) N=N+1 
10 CONTINUE 
II=1 
Bl IF(PCPN(5).EQ.R.OR.PCPN(5).EQ.C) GO TO 23 
JJ=0 
24 DOW J=1 34 
IF(PCPN(J) .EQ.R.OR.PCPN(J) .EQ.C.OR.PCPN(J) . EQ. BLANK) 
1GOP TO) 15 
i CONTINUE 
DO 12 M=1,5 
IF(PCPN(M).EQ.T) PCPN(M)=ZERO 
12 CONTINUE 
21 DO 70 K=1,5 
WRITE (77102) ID(K),1Y(K),IMCK)) (T1(K,1) ,1=1,9), PGPN{K) , 
LCT2 (hs Le i=1 68) 
70 CONTINUE 
41 IF(JK.EQ.3) GO TO 28 
II=II+1 
IF(I1.EQ.6) GO TO 18 
DO 19 K=1,5 
19 READ(4,102) ID(K),1Y(K),IM(K), (T1(K,I) ,1=1,9),PCPN(K), 
LGR ke 1 je isl,8) 
LE(JISEQe1) .GO=TO 20 
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APPENDIX D (CONT'D) 


CORTON 1 

JJ=1 

DO 16 J=1,5 

PCPN(J) =DNEG 

COLTON? 

IF(PCPN(5).EQ.R.OR.PCPN(5).EQ.C) GO TO 25 

JJ=0 

DO 22 J=1,5 

PCPN(J)=DNEG 

GOrT0 +21 

JJ=1 

GO TO 26 

IQ=N-27 

GO TO (1,2,3,4),1Q 

DO 27 K=1,3 

READ(4,102) ID(K),IY(K),IM(K), (T1(K,I) ,I=1,9) ,PCPN(K), 
1(T2(K,1I),I=1, 8) 

GO TO 28 

DO 29 K=1,4 

READ(4,102) ID(K) ,IY(K),IM(K) , (T1(K,1I),1=1,9) ,PCPN(K), 
1 (P2¢KS1)) 1=1,8) 

GONrTO 28 

DO 30 K=1,5 

READ(4,102) ID(K),IY(K) , IM(K), (T1(K,I) ,I=1,9),PCPN(K), 
1(T2(K,1) , T=1,8) 

JIK=3 

G0 TO -31 

DO 32 K=1,6 

READ(4,102) ID(K),IY(K) , IM(K) , (T1(K,I) ,I=1,9) ,PCPN(K), 
1(T2(K,I) , 1=1, 8) 

JIK=3 

GO TO 31 

STOP 

END 
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APPENDIX E 


Program for Calculating Characteristic Curve Parameters for 
Ordinary Climatological Stations 


FORTRAN 


DIMENSION A(31,12,6) ,B(31,12,6),PCPN(5) ,NAME(6), 
1LMDA1(6) , LMDA2 (6) , IY(5) , IM(S) 
REAL LMDA1, LMDA2 

100  FORMAT('4',45X, "PROGRAM FOR ORDINARY CLIMATE STATION'//) 

101  FORMAT(I7,6X,6A4) 

102 FORMAT(7X,212,33X,F4.2) 

103 FORMAT(50X,17,4X,6A4,//) 

104 FORMAT(14,6(2X, F6.3,2X,F6.3)) 

105  FORMAT('MNTH',2X,'LMDA1', 2X, 'LMDA2' ,4X, 'LMDA1', 2X, 'LMDA2', 
14X,'LMDA1',2X,'LMDA2',4X, 'IMDA1', 2X, 'LMDA2', 4X, 'LMDA1', 2X, 
2"LMDA2', 4X, 'LMDA1', 2X, 'LMDA2'/) 

WRITE (6, 100) 
READ (7,101) IP, (NAME(I),I=1,6) 
WRITE (6,103) IP, (NAME(I),I=1,6) 
WRITE (6,105) 
IFY=26 
IFM=1 
II=1 
IPY=1 
Z=0.0 

4 DO 14 K=1,5 
READ (7,102,END=8) IY(K),IM(K) , PCPN(K) 
IF (PCPN(K).GT.0.0) Z=Z+1 

14 + CONTINUE 

34 IF (IM(1).NE.IFM) GO TO 9 
IF(IY(1).NE.IFY) GO TO 36 
SUMA=0.0 
DO 5 K=1,5 
SUMA=SUMA+PCPN(K) 

5 CONTINUE 

A(IPY, IFM,I1)=SUMA 

B(IPY, IFM,II)=Z 
II=II+l 
Z=0.0 
IF(IFM.EQ.2.AND.II.EQ.6) GO TO 19 
IF(II.EQ.7) GO TO 19 
= GOmtO 4 

19 IFM=IFM+1 
II=1 
Z=0.0 
IF(IFM.GT.12) GO TO 16 
GO TO 4 

16 IFY=IFY+1 


36 
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18 


ool 
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APPENDIX E (CONT'D) 


IPY=IPY+1 
IFM=1 

II=1 

Z=020 

(Pelee hoje GOsd als 
GO TO 4 

DO 13 IFM=1,12 

DOpe oe 21s 6 

A(IPY, IFM,II)=-.99 
BULEN aL EM I= =. 99 
CONTINUE 

CONTINUE 

IFY=IFY+1 

IPY=IPY+1 

IFM=1 

Ii=1 

7=020 

IF(IFY.GT.56) GO TO 18 
GO TO 34 

NOnoo ail 6 
NCIPY .LEMe LL) =—.99 
B(IPY, IFM, II)=-.99 
CONTINUE 

IFM=IFM+1 

If=1 

720.0 

IF(IFM.GT.12) GO TO 44 
GO TO 34 

One LEM= Ll 
IF(IFM.NE.2) IMAX=6 

IF (IFM.EQ.2) IMAX=5 

DOs) 7a) l=t, IMAX 
COUNT=0.0 

SUML1=0.0 

SUML2=0.0 

De te LY = les 

TPCAT LEY 1FMaIL). LT. 0.0)"GOe10 21 
SUML1=SUML1+A (IPY, IFM, II) 
SUML2=SUML2+B(IPY, IFM,II) 
COUNT=COUNT+1 . 
CONTINUE 

LMDAI (I1)=SUML2/ (5.0*COUNT) 
LMDA2 (II) =SUML2/SUML1 
CONTINUE 

WRITE(6,104) IFM, ((LMDA1(II), LMDA2(II)),1I=1, IMAX) 
CONTINUE 

STOP 

END 
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APPENDIX F 


Program for Solution of Distribution Function of X 


Ca aa eS eee eran | 


FORTRAN 


C PROGRAM XSUBT 
DIMENSION T(6),FF(100) , FFF (100) 
LAMDAS=1 
DX=0.05 
READ(5,4) T 
4 FORMAT (6F10.1) 
DO 50 M=1, LAMDAS 
READ(5,5) AL1,AL2 
5 FORMAT (2F10. 3) 
DO 50 K=1,6 
FTX2=0.0 
WRITE(6,51) AL1,AL2,T(K) 
51) FORMAT (IH1,'LAMDA1 = ",F5.3,'LAMDA2 = ',F5.2,' T = ',F5.1) 
WRITE(6,52) 
52. FORMAT(1H,' X DIST. FUNT. INCREMENTS' ) 
X=0.0 
A=AL1*T(K) 
FTX 1=EXP (-A) 
SUMX=FTX1 
GOwr0"54 
3 B=AL2*X 
SUML=0. 0 
DO 60 J=1,50 
AJ=J 
SUML=SUML+ALOG (AJ) 
FFL=-B+J*ALOG (B) -SUML 
FE (J) =EXP (FFL) 
FFF L=-A+J*ALOG (A) -SUML 
60 FFF(J)=EXP(FFFL) 
SUM1=EXP (-B) 
DO 12 J=1,50 
12. SUM1=SUM1+FF (J) 
SUMX=SUM1*EXP (-A) 
DO 15 L=1,50 
SUM=0.0 
DO 14 J=L,50 
14. SUM=SUM+FF (J) 
SUM=SUM*FFF (L) 
15. SUMX=SUMX+SUM 
FTX1=SUMX 
IF (ABS(1.0-FTX1) .LT.0.001) GO TO 50 
54.  DETX=FTX1-FTX2 
IF (X.EQ.0.05) DFTX=FTX1 
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APPENDX F (CONT'D) 


WRITE (Oso) eee TX 1 DELX 

FORMAT (SX ,F5.2,6X,F9,3;6X,F9.3) 
FTX2=FTX1 

X=X+DX 

CONTINUE 

STOP 

END 





Meteorological Information of the Selected Stations 


Station Name 
Cardston 
Carway 
Coleman 
Edmonton 
Edson 
Entrance 
Fort MacLeod 
Jasper 

_ Lethbridge 
Lundbreck 
Medicine Hat 
Seven Persons 


Winnifred 


beh a. ak! 
Ae aed 
49°00! 
49°38! 
seh sel 
53°35! 
CS 28! 
49°43! 
52°53! 
49°38! 
AO255)! 
50.01! 
49°50! 


49°54! 
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APPENDIX G 


TABLE 1 


Long. W 
Lise tos 
psa 
IAS S5! 
GS Sal" 
Tie, 25! 
11742) 
115524! 
118°04! 
11a Ace 
114-08! 
110243! 
10.54 


Tel 2 


+ Present elevation of the station 


Elevationt 


SUAS 


4460 


4400 


poe 


S0Ss 


3300 


3116 


3480 


3018 


3918 


Zao) 


2480 


LT ao 


++ Period of record chosen for this study 


SELECTED PRECIPITATION STATIONS OF ALBERTA 


Years of 


Recordtt 


1956 


1956 


1956 


1968 
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APPENDIX G (CONT'D) 


iyi 


TABLE II 


PRECIPITATION DATA FOR THE SELECTED STATIONS OF ALBERTA 


tation Annual 
Name Precipitation 
Cardston 18.04 
Carway PAO ASS 
Coleman Oe L 
Edmonton 18.64 
Edson 20.85 
Peconic e 1o.96 
Fort MacLeod 17.67 
Jasper 15296 
Lethbridge 17.25 
Lundbreck 19.28 
Medicine Hat 14.29 
Seven Persons 13.79 
Winnifred 135286 


Note: The vaiues 


presented in this table are based on the 


Annual 


Rainfall 
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ike 


Toe 


14, 


14, 


IPAe 


i: 


LO. 


BOs 


8. 


ROU 


&6 


21 


26 


98 


4] 


38 


. 06 


66 


OS 


Annual 
Snowfall 


65. 


OD 


Cie: 


ee 


D0% 


Don 


Dae 


AQ, 


G3. 


San 


48. 


A7. 


49, 
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0 


6 


4 


No. of Days 
With Rain 


47 
43 
46 
68 
69 
o4 
oo 


§1 


45 


period of observation from 19350 to 1960 


No. of Days 
With Snow 


47 
43 
34 
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APPENDIX H 


Computational Method for the Lambda Calculations 


TABLE@ iid 
COMPUTATIONAL METHOD USED TO DERIVE THE AX VALUES 
Station Carway, Alberta 
Month January 
Number of Storms 


Period (Days) 


Year 1-5 6 - 10 11 - 15 16 - 20 Zi 25 26 - 30 
1926 0 3 0 etc 
1927 2 2 etc 
! ' 1 
1956 1 etc 
Total 28 etec 


E(n,) = 28/29 = .966 etc 


r 





= = ‘SO = 
l rs .966/5 LOS etc 


where 29 is equal to the number of years of record and t is the interval 


length. 
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APPENDIX H (CONT'D) 


TABLE SLY, 
COMPUTATIONAL METHOD USED TO DERIVE THE A, VALUES 
Station Carway, Alberta 
Month January 
Total Amount of Precipitation 


Period (Days) 


Year 1-5 6 - 10 ise ES 16 - 20 2h 25 26 - 30 
1926 0 co2 0 etc 
1927 200) 226 Ee. 
' ' ' 
1956 . 06 etc 
Total 4.85 etc 


E(X,) SEAS 57 290s 8. G67 etc 
re eS 87s etc 


where 29 is equal to the number of years of record and t is the interval 


length. 
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APPENDIX I 


Comparisons Between the Observed and Theoretical Distributions 
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TABLE VII 


EDMONTON - NO. OF DAYS WITH PRECIPITATION DURING 5-DAY PERIODS 





No. of days 
with 

precipitation 0 i 2 3 4 5 
Jaleo I = 0.388 

Obsvd. freq. 8 6 10 3 3 i 
iieor erred. 5 S) 8 > 35 1 
Feb, 1-5 Ay = 0.383 

Obsvdes freq: 5 4 12 6 4 0 
Theor y-freq: 5 9 8 5 3 1 
Mar. 1-5 Ay = 0.360 

Qbsvd. freq. = 8 is) 5 1 1 
wheor. freq. 5 9 9 i) 2 ai 
fpr 1-5 AY = 0.275 

Obsvd. freq. 10 6 6 8 1 0 
Theor, freq. 8 if 1 4 q. 0 
May 1-5 Ay = 0.258 

Obsvd. freq. 8 8 i) = 1 1 
TROT ue GCC. 9 val vf 3 1 0 
June 1-5 Ny = 0.335 

Obsvd. freq. 4 10 WH 3 Z 0 
Theor. freq. 6 10 8 4 2 1 
July 1-5 Ay = 0.440 

Obsvd. freq. 3 8 8 8 4 0 
ff eCOireEGed . 3 8 8 6 4 Z 
Aug. 1-5 Ay = 03595 

ODSVdemt req. 4 5 8 13 i 0 
WReor er Led. 4 3 8 6 5 1 
Sept. 1-5 hy =40 5552 

Obsvd. freq. 5 15 4 4 3 0 
Tico, tred> 6 10 8 4 2 1 
Cct. 1-5 Ay = 0.246 

Obsvd. freq. 8 8 a 1 2 0 
TiCOc se ire. 9 11 v 3 1 0 
Nov. 1-5 Ay =e e2as 

Obsvd. freq. MS 9 6 2 1 0 
Theoteg ire. 10 il if 2 1 0 
DEC. iso AY =s 02515 

Obsvd. freq. 8 7 12 2 2 0 


TNeOnem LEC 6 10 8 4 4 1 
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TABLE VIII 


CENTRAL ALBERTA STATIONS - NO. OF DAYS WITH PRECIPITATION 
DURING 5-DAY INTERVAL JAN. 1-5 


No. of days 
with 
precipitation 0 1 2 3 S 5 
a erie ee er ct en i ea ae et a gh a ie nce 


Edmonton AY = 0.388 


Obsvd. freq. 8 6 10 3 2 val 
Pheorasrred. 5 9 8 5 3 1 
Edson Ay = 0.212 

Obsvd. freq. 16 Z 5 2 i 0 
PHheO rast red. it Tt 6 2 i 0 
Entrance AY ca Oe 

Obsvd. sireq: 14 6 5 3 2 0 
Theor seired. 10 11 6 Z 1 0 
Jasper Ay = 0.350 

Obsvd. freq. 4 9 6 8 2 1 
THeOr ee cred: 5 9 8 5 2 1 
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TABLE IX 


EDMONTON AND EDSON - NO. OF DAYS WITH PRECIPITATION - PERIOD JUNE 1-20 





No. of days Edmonton Edson 
with Ne = 0.961 A= 0.421 
precipitation Obsva-= freq. Iheor. freq. Dbsvdired . aheor: freq. 
0 0 0 0 0 
il 0 0 0 0 
2 i 0 0 0 
3 2 1 Z st 
4 1 y Z 2 
5 3 3 i 2 
6 3 4 5 3 
fe 4 5 2 4 
8 i 4 » 4 
9 4 4 2 4 
10 4 % 1 = 
11 4 Z 5 3 
12 4 1 3 2 
Ls 0 1 2 1 
14-20 0 1 1 Z 
TABLE X 


ENTRANCE AND JASPER - NO. OF DAYS WITH PRECIPITATION - 
PERIOD MAY 21 - JUNE 10 





No. of days Entrance Jasper 
with hee US 54 A, = 0.349 
precipitation Obsvd.gireq.. [heor mereq. Obsvd wired. aiheor.gtredq. 
0 0 0 0 0 
1 0 0 0 0 
2 i i 0 1 
3 2 1 2 2 
4 3 3 2 3 
5 J 5 7 4 
6 5 4 3 5 
7 3 4 4 5 
8 4 4 5 4 
9 2 3 3 3 
10 Ss Z 2 2 
il i 1 4 1 
12 1 1 0 1 
15 1 f 0 0 
14-20 1 0 1 0 
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TABLE XI 


JASPER AND EDSON - NO. OF DAYS WITH PRECIPITATION - PERIOD JAN. 1-30 





No. of days Jasper Edson 
with Ay = 07350 Ay =s05212 
precipitation Obs Vee ered a heer. sane. Obsvd.) freq. Theor. sfreq;, 
0-1 0 0 1 0 
2-3 1 0 8 3 
4-5 5 1 3 8 
6-7 5 4 6 10 
8-9 5 6 6 6 
10-11 3 i 6 3 
12-15 4 6 0 1 
14-15 5 5 il 0 
16-17 5 2 0 0 
18-19 1 1 0 0 
20-30 0 0 0 0 
TABLE XII 
EDMONTON AND ENTRANCE - NO. OF DAYS WITH PRECIPITATION - 
PERIOD JAN. 1-30 
No. of days Edmonton Entrance 
with hy = 0.385 AD =a Oeil 
precipitation Obsvd. -£red, ml heor sired. ODsVd. UGedw NeOr.erred: 

0-1 0 0 Z 0 
2-3 1 0 8 5 
5 2, 1 6 8 
6-7 3 3 2 9 
8-9 4 5 6 6 
lo-11 6 7 4 5 
12-13 5 7 0 1 
14-15 4 5 0 0 
16-17 4 2 0 0 
18-19 2 i 2 0 
20-30 0 0 0 0 
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TABLE XIII 
CENTRAL ALBERTA STATIONS - AMOUNT OF PRECIPITATION DURING 5 DAYS 
te =O NnUuSs. wane eS ce 
frequencies 


Amount of Edmonton Edson Entrance Jasper 
precipitation Obsvd. Theor. Obsvd. Theor. Obsvd. Theor. Obsvd. Theor. 





0 8 4 16 gs! 14 10 S S 

0 = 0720 12 Ly 8 id 11 WZ 16 13 
0.21 - 0.40 9 iy 2 5 5 5 3) 7 
0.41 - 0.60 i 2 2 2 0 2 0 3 
0.61 - 0.80 1 1 I 1 0 1 2 it 
0.8 25-2200 0 0 0 1 1 0 2 a 
Ol = E.20 0 0 2 0 0 0 0 0 
Poal, — 40 0 0 0 0 1 0 0 0 
1.41 - ©& 0 0 0 0 0 0 1 0 

TABLE XIV 


CENTRAL ALBERTA STATIONS - AMOUNT OF PRECIPITATION DURING 10 DAYS 
bs = 0 hrs. Jairpeerse 
frequencies 


Amount of Edmonton Edson Entrance Jasper 
precipitation Obsvd. Theor. Obsvd. Theor. Obsvd. Theor. Obsvd. Theor. 


0 1 1 6 4 Ht 4 it 1 

0. = 0.30 18 16 16 14 Ls 16 Ls Lz 
Cale — 0560 10 10 6 8 2 % 8 10 
0.61 - 0.90 1 3 1 4 a 2 1 J 
Oa Se = eZ0 0 1 1 1 1 1 4 Z 
Bowe — Yas0 | 0 1 0 Zz 0 1 0 
Gaok = 8480 0 0 0 0 0 0 i 0 
ES 2710 0 0 0 0 0 0 0 0 
sac ed Soueraial d 0 0 0 0 0 0 24 0 
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TABLE XV 
CENTRAL ALBERTA STATIONS - AMOUNT OF PRECIPITATION DURING 15 DAYS 
a == Or nus.. Jan, lst 
frequencies 


Amount of Edmonton Edson Entrance Jasper 
precipitation Obsvd. Theor. Obsvd. Theor. Obsvd. Theor. Obsvd. Theor. 


0 0 3 1 5 1 0 0 

0 - 0.40 16 14 14 14 iF, 16 2 10 
0.41 - 0.80 11 13 ~) 10 4 9 a ih 
Ooi 20 3 =) 3 4 i! 3 fi 6 
¥.20 -— Y..60 1 il 2 2 I 1 3 Z 
Gk = 4.00 0 0 0 0 1 0 I 0 
2.01 - 2.40 0 0 0 0 i 0 0 0 
At — 2.80 0 0 0 0 0 0 0 0 
2s <3 0 0 0 0 0 0 0 1 0 

TABLE XVI 


CENTRAL ALBERTA STATIONS - AMOUNT OF PRECIPITATION DURING 20 DAYS 
te =F OF bse, Jans ise 
frequencies 


- Amount of Edmonton Edson Entrance Jasper 
precipitation) Obsvdi. theor., Obsvd. Theor. Obsvd., Theor. Obsvd.. Theor. 


0 0 0 0 0 1 0 0 0 

0 - 0.40 11 8 is 10 15 ie 8 5 
0.41 - 0.80 14%) iS = Tt 7 i] 6 a 
O.. 8) =. 1.20 6 6 8 6 4 5 6 a 
21h 1 -12..G0 0 2 2 3 ‘| i 8 4 
Git =. 27.00 1 0 1 ih 1 1 1 1 
2.01 -. 2.40 0 0 0 0 1 0 0 0 
2.41 - 2.80 0 0 0 0 0 0 0 0 
Dee 36 20) 0 0 0 0 0 0 0 0 
Se 2A, ns 0 0 0 0 0 O- 0 0 
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TABLE XVII 
CENTRAL ALBERTA STATIONS - AMOUNT OF PRECIPITATION DURING 25 DAYS 
Ce = Oshrs.. Alan. absits. 
frequencies 


Amount of Edmonton Edson Entrance Jasper 
precipitation, Obsvdky Theor. Obsvd,. Theor. Obsvdee Theor). Obsvd.e Theor; 


0 0 0 0 0 1 0 0 0 

0 - 0.50 Te x 14 2) 11 12 9 4 
Deo ALs00 sO, 15 6 Me 10 12 8 1 
PSO 4= 41, 450 8 6 6 iz, 4 5 9 2 
Lito s 42/00 1 2 4 2 2 1 S 4 
Zeal = #2;.500 0 0 0 1 2 0 i 1 
Pato =) 5.00 0 0 1 0 0 0 0 0 
STONES SF 1) 0 0 0 0 0 0 0 0 
Bwbl = 00 0 0 0 0 0 0 1 0 


TABLE XVIII 
CENTRAL ALBERTA STATIONS - AMOUNT OF PRECIPITATION DURING 30 DAYS 
Lae Ogre. Jan. wl Site. 
frequencies 


Amount of Edmonton Edson Entrance Jasper 
precipitation Obsvd. Theor. Obsvd. Theor. Obsvd. Theor. Obsvd. Theor. 
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TABLE XIX 
MEDICINE HAT AND SEVEN PERSONS - AMOUNT OF PRECIPITATION DURING 5 DAYS 
Ue = 0 hrs Aug. lst. 
frequencies 


Amount of Medicine Hat Seven Persons 
preerprcacvon » Obsvdes freq. Theor. fred. Obsvda fred. Ineor.- freq. 


0 13 g 19 15 

OF = Omz0 14 10 4 6 
Omar = OF40 1 6 4 4 
0.41 - 0.60 3 0 2 
0.61 - 0.80 0 Zz 2 2 
Oacie— Ies00 0 1 1 i 
BeOls = 20 1 0 a 1 
1.21 - 1.40 i 0 0 0 
1.41 - @ 0 0 0 0 

TABLE XX 


MEDICINE HAT AND SEVEN PERSONS - AMOUNT OF PRECIPITATION DURING 10 DAYS 


Ue = O'hrsy Augie Ist. 


frequencies 


“Amount of Medicine Hat Seven Persons 
preciprtaciom  Obsvdy freq Theor. freq). Ubsvd. fred. Theor. fred. 





0 4 3 10 a 

0 - 0.30 Th tee 8 10 
Os) O60 4 i, & 6 
OsGlr =" 0790 th 5 = 4 
GaG Pe 120 1 5 2 Z 
Pea = 50 0 i 1 a 
ies = 1280 2 1 Z 1 
eo = 200 2 0 0 0 
ogg 8 ai a 0 0 0 0 
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TABLE XXI 
MEDICINE HAT AND SEVEN PERSONS - AMOUNT OF PRECIPITATION DURING 15 DAYS 
te a0 nes. EAs sc. 
frequencies 


Amount of Medicine Hat Seven Persons 
Drees DAtati on | Obed. req, (Theor. freq. | \Obsyvd. atreq.. Theor. freq, 


0 3 1 a 2 

0 - 0.40 14 10 8 1l 
0.41 - 0.80 5 10 5 8 
Uso) =) 1.20 “ 6 4 5 
i221 1.00 2 & 4 2 
1GOl = 2.00 i 1 % | 
2.01 - 2.40 a 0 0 1 
241 - 2.80 i 0 0 0 
2.01 - ih 0 1 0 

TABLE XXII 


MEDICINE HAT AND SEVEN PERSONS - AMOUNT OF PRECIPITATION DURING 20 DAYS 


Ce = 90 vba s | GAUC asic 





frequencies 
~ Amount of Medicine Hat Seven Persons 
misecipacation iObsvd.“freq. Theor. freq. Obsvd, sireq. Gheor.yoneq: 
1 0 6 2 
0 - 0.40 tel 6 tf 8 
0.41 - 0.80 6 9 3 8 
O81 - 1.20 4 8 5 6 
ia = (1.60 4 5 4 5 
AotGn 222100 1 2 5 2 
2201 - 2.40 Z 1 0 1 
2441 .- 2.80 0 0 0 1 
PE, Re PAM) 1 0 0 0 
BdiZdun © il 0 i 0 
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TABLE XXIII 
MEDICINE HAT AND SEVEN PERSONS - AMOUNT OF PRECIPITATION DURING 25 DAYS 
ue =90) hrs'rtAugiel st. 


frequencies 


Amount of Medicine Hat Seven Persons 
precipitation Obsvd. freq. ‘Theor. freq. Obsvd setreq.. Theor, trea. 





0 1 0 4 1 

0 = 0750 10 5 6 8 
0.51 - 1.00 5 10 6 ) 
Ol = 50 5 8 2 6 
Peak = °27.100 2 5 8 4 
rie =" Aa 0 4 2 & Z 
Se 3 OO 2 1 1 1 
ade <P StU 1 0 0 0 
Set be = oe iL 0 I 0 

TABLE XXIV 


MEDICINE HAT AND SEVEN PERSONS - AMOUNT OF PRECIPITATION DURING 30 DAYS 


- =e ONS taeAUg seal Site 


frequencies 
“Amount of Medicine Hat Seven Persons 
Puceipiracione  Obsvd. freq. i!heor. freq. Obsvd™etreq. -Theor, tired. 

0 0 3 0 
OF —* 07,90 10 3 6 6 
Oho ~— 14.60 5 8 4 8 
PO — 150 5 ) 2 rl 
DSi, = 2.00 3 6 8 s) 
Z0l, = 2.50 a 3 4 5 
aisle = Svs 5 1 0 1 
Se = 3.50 u 1 2 1 
3,51) -- 4700 0 0 2 0 
da} Leecoereico 2 0 0 0 
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TABLE XXV 
AMOUNT OF PRECIPITATION DURING 5 DAYS FOR CENTRAL ALBERTA 
Average frequencies of the central Alberta stations combined 


Le SAG) jnteay Aine dies 


Amount of Obsvd. Theor. Theor. 
precipitation ave. ave.) ave. , 
0 10 8 1. 

0 - 0.20 ie 13 14 
O-241 = 0740 5 6 6 
0.41 - 0.60 il 2 2 
ete Sete) 1 1 1 
Oe Osa 00 il 1 1 
VeOL —- 61.20 1 0 0 
Le2k i 40 0 0 0 
Tl. 4) =)" (0 0 0 0 

TABLE XXVI 


AMOUNT OF PRECIPITATION DURING 10 DAYS FOR CENTRAL ALBERTA 
Average frequencies of the central Alberta stations combined 


cs =e Mri Shercll. am lS te 





“Amount of Obsvd. Theor. Theor. 

precipitation ave ave.) ave. 
0 5 3 2 
0 - 0.30 NBS 15 15 
0.31 - 0.60 7 9 9 
0.61 - 0.90 1 5 4 
U,9l = 1.20 2 1 1 
Leelee LU 1 0 0 
[51 - 1.80 0 0 0 
Poke = 2k) 0 0 0 
Zs O01 <= 0 0 0 
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TABLE XXVII 
AMOUNT OF PRECIPITATION DURING 15 DAYS FOR CENTRAL ALBERTA 
Average frequencies of the central Alberta stations combined 


a =A enwS,modn. sists 





Amount of Obsvd. THeor, Theor. 
precipitation ave, ave.) ave 
@ it 0 
0 = 0.40 15 14 14 
0.41 - 0.80 7 ia 1 
Oe81- = shs20 4 4 4 
lara = alk 60 Z 1 1 
i641. == 27,00 1 0 0 
2.01 — 24/40 0 0 0 
2041 = +2,80 0 0 0 
Zeol = 0 0 0 


TABLE XXVIII 
AMOUNT OF PRECIPITATION DURING 20 DAYS FOR CENTRAL ALBERTA 


Average frequencies of the central Alberta stations combined 


Ue = "0 shrs, Jan. ast. 
Amount of Obsvd. Theor Theor. 
precipitation ave. ave., ave. 
0 0 0 0 
0 - 0.40 12 9 8 
0.41 - 0.80 8 Le 13 
0.81 = 1.20 6 6 7 
deel = 15,60 3 2 2 
1.61 = 2.00 1 1 1 
2,01 = 2.40 0 0 0 
2.41 - 2.80 0 0 0 
ool ee Se 20) 0 0 0 
Se ba 0 0 0 













t. WY RATTATTS VENt 3 — . "— 


+] 
—_— | tk tines St Fo estomenge signa 
; a & 


Wi 


7 a - 
f ‘ % — _ 
° “ te 
{ : j ff 


: as 
, baiy - to anvomA — 
‘\ TER 


= ae - _ 

| a1 a 
ee ; -_ = O80 Theo 

; - ‘osct. = 18.9 7 

: Ga. i = AS 





=. 


A 
“wis to 
Gh .5 
AS 
o- 
04.4. =< I 
i a . 
; 
‘ 
2 
> oy : 7 J " 
| PPATP + FO TAUOMA 
& © 
\ ) 7 
: ‘3 “eens Tt. SyBrTOVA 
oa ae. 7 
- 


1) i : | ae ee 


He .0 a 


15 


TABLE XXIX 
AMOUNT OF PRECIPITATION DURING 25 DAYS FOR CENTRAL ALBERTA 
Average frequencies of the central Alberta stations combined 


ue =SmOenirs moan USt. 





Amount of Obsvd., Micon. Theor. 

precipitation ave. ave., ave. 
0 0 0 0 
OR =e 07 50 at 8 7 
Qeol aL. 00 5 13 14 
i lB ee NO) Zc Ul 7 
Ve ees 2 00 2 2 2 
AA” PABEAD if 1 0 
2 ole 00 0 0 0 
SWS Seb 0 0 0 
eo em aco 0 0 0 

TABLE XXX 


AMOUNT OF PRECIPITATION DURING 30 DAYS FOR CENTRAL ALBERTA 
Average frequencies of the central Alberta stations combined 


ty =sQenrs.s. Jan. LSt. 


Amount of Obsvd. Theor. Theor: 

precipitation ave. ave.) ave. , 
0 0 0 0 
OF- ae 50 7 S 4 
Os eral, OO ital i2 a) 
OL — Loo0 6 8 9 
lee OU 5 4 4 
2.01 - 2.50 ih L 1 
PSN IUD 1 1 0 
BG AO IE Bi ENO 0 0 0 
Se les 00 0 0 0 
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Fig 6. , values computed using distribution of n, 
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Fig 8. , values computed using distribution of ne 
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Fig.10. d, values computed using distribution of My 
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Fig.il. X, values computed using distribution of ny 
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PERIOD 1926 — 1956 
COLEMAN, ALBERTA 


Time of Year 


Fig.!2. d, values computed using distribution of n 
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Fig.13. », values computed using distribution of n 
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Fig.19. X 2 wolues computed using distribution of X, and Xy (vt) 
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Fig.22. 2g volues computed using distribution of X, and X.(v=l) 
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Fig.24. 2 values computed using distribution of X, and. Xy(v=l) 
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Fig.25. 2, values computed using distribution of X, and Xy (vel) 


22 mise ' esnin “ 


Pn ae hd i eade 
hall ho ae is ; 4 A 7 — D i ms - 


er 5. B.. = 


os 





= ne he pte nf 


on A conan ( : om >... i 
/ | iL > : © ? . - : i} ‘ oO ' 7 ; 
' 





| : | eh 
; . } “ , . es) 
= Gq ; : : ; : | . LF 
; = ee oe ee ee ~ - b—- — i | —_+ - — = ok gine EE a 
| + } ; ; } ; ' = le. 
Pb / iis . . at 
s : | 
ape i } } | J a ' - : | - 
>) if 4 } . 
_ ee _ pacd — = "| ; : a a 
_— AD 7 ‘1 
7 7 | i | ] 6 
; ; 


d2 Values (I\/amount of pcpn per storm) 


2 Values (i/amount of pcpn per storm) 


PERIOD (938 - 1968 


EDMONTON, ALBERTA 


Choracteristic 2 curve 


| ooo 
see IES EIN 





Time of Year 
Fig. 26. 2 values computed using distribution of X, and Xytve 


PERIOUMIS 23959 
EDSON, ALBERTA 


Characteristic X, curve 





JAN FEB MAR APR MAY JUN JUL AUG SEPT OcT NOV DEC 


Time of Year 


Fig. 27. », values computed using distribution of X; and Xy (v=1) 
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Fig. 36. Comparison of the characteristic X, curves for 
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Fig.38. Characteristic X, curve with first standard deviation for 
the southern Alberta stations 
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Fig.40. Comparison of the characteristic X, curves for 
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Fig.42. Characteristic XX, curve with first standard deviation 
for the central Alberta stations 


Characteristic %» Curve 
Central Alberta Are emmrama= 


Southern Alberta Area ous am 


MAR APR MAY JUN JUL AUG SEPT OcT NOV DEC 
Time of Yeor 
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Fig,46. Comparison of the characteristic », curves for 
Jasper and Coleman 


Characteristic > curves 
JOSDCl cece sa 


Coleman aus ous ane axe axe axe ae 


MAR APR MAY JUN JUL AUG SEPT OcT NOV DEC 
Time of Yeor 


Fig.47. Comparison of the characteristic >, curves for 
Jasper and Coleman 
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Fig. 48. Comparison of the characteristic 2», curves for 
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Fig. 49. Comparison of the characteristic >, Curves for 
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Fig. GO Edmonton~ Average monthly no. of days with measureable 
precipitation 
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Fig.63 Seven Person - Average monthly amount of precipitation 
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